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FINAL  SCIENTIFIC  REPORT 

TITLE:  THEORETICAL  STUDIES  OF  ATOM-SURFACE  INTERACTIONS 

Grant  No.  ONR  N00014-75-C-0949 

Principal  Investigator:  Professor  Narkis  Tzoar 

The  City  College  of  The  City  University  of  New  York 

The  main  thrust  of  our  research  was  to  investigate  the 
desorption  of  atoms  from  solid  state  surfaces.  Interest  in  atom 
desorption  stems  from  several  sources.  First  it  provides  a  tool  for 
characterizing  surfaces.  One  may  identify  and  measure  the  amounts 
of  various  adsorbates  on  a  surface  and  learn  the  elementary  physical 
and  chemical  processes  associated  with  their  bonding.  Secondly,  as 
a  technological  tool,  desorption  provides  technique  for  preparing 
and  cleaning  a  surface.  Thirdly,  the  degradation  of  the  surface  and 
underlying  substrate  of  the  walls  by  radiation  or  particle 
bombardment  become  a  matter  of  practical  concern  in  many  situations. 
Also  desorption  of  atoms  from  the  surfaces  of  interstellar  dust 
grains  provides  a  mechanism  for  the  coupling  of  radiant  energy  from 
stars  to  the  kinetic  energy  interstellar  atoms  or  molecules. 
Desorption  is  thus  a  phenomenon  of  rather  universal  concern  and  of 
practical  interest  in  many  branches  of  physics  and  chemistry. 

Experimental  and  theoretical  developments  related  to  this 
project  were  reviewed  recently,  in  the  November  1980  issue  of 
Physics  Today  (see  Search  and  Discovery  column). 

One  may  regard  desorption  as  a  collision  process  in  which  an 
incident  projectile  collides  with  an  adatom  and  ejects  it.  If  the 
incident  particle  is  an  ion  the  process  is  known  as  sputtering.  If 
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the  incident  particle  is  an  electron  or  photon,  it  termed  electron 
stimulated  desorption  or  photodesorption  respectively. 
Thermaldesorption  may  be  thought  oE  as  a  collision  process  involving 
incident  thermally  excited  phonons. 

The  main  efforts  during  the  first  half  of  the  grant  duraction 
was  two  fold:  (a)  to  develop  the  kinetic  theory  of  the  ion 
distribution  in  electron  stimulated  desorption  (ESD) ,  (b)  to 
determine  the  potential  energy  curves  which  repel  the  ions  away  from 
the  host  surface. 

Our  kinetic  theory  assumes,  for  simplicity,  a  two  level  system 
i.e.  an  energy  curve  for  the  adatom  which  determines  its  ground 
state  and  an  excited  (repulsive)  energy  state  for  the  adatom's  ion. 
The  parameters  which  determine  these  energy  curves  will  be 
considered  shortly.  We  then  obtain  the  coupled  Boltzmann  equations 
for  the  deviation  of  the  probability  from  equilibrium  for  finding 
the  adatom  in  one  of  the  two  states  as: 

,P  •  (1) 

^  +  v.pp, -1-t ?er§ 

Here  E^  (E2)  is  the  energy  curve  of  the  ground  (excited)  state 
and  the  subscripts  "E"  and  "Q"  stands  for  "Excitation"  (by  the 
incoming  electron)  and  "Quenching"'.  The  source  term  (P2)E 
taken  to  be  weak  thus: 

0 i)e*  Po* 

Where  R  is  a  rate  given  from  scattering  theory  and  PQ  represents 
the  ground  state  of  the  adatom.  The  quenching  term  is  taken  to  be 


W  )  tfl'v 
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(3) 


Here  Tq  is  an  average  quenching  time  on  the  substrate  surface 

and  ?  represent  the  range  of  the  "substrate"  electrons  away  from  the 

surface.  We  thus  may  write  the  kinetic  equation  as: 


Hi.  +v.vp 

3-t  2 

+  v  .  V  p, 

7>t  ' 


(4) 


This  set  of  equations  yields  the  angular  and  energy  dependence, 
linear  in  the  incident  flux,  in  terms  of  the  energy  curves  (E^ , 

E2) ,  the  excitation  rate  (R)  and  the  quenching  (Q) . 

In  view  of  the  complexities  associated  with  the  solution  of  Eq.  (4) , 
we  adopt  an  approach  well  suited  to  numerical  simulation  of  ESD  ion 
angular  and  energy  distributions. 

The  steady  state  ESD  response  is  found  by  solving  an  equivalent 
time  dependent  problem.  We  simulate  the  flux  of  ions  or  neutral 
atoms  desorbed  in  response  to  an  incident  electron  pulse  of  short 
duration  X  ,  then  take  the  limit  T. — >  0.  We  obtain  below  suitable 
expressions  for  the  total  desorption  cross  sections,  using  the 
transport  formulation.  Numerical  differentiation  in  the  asymptotic 
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plane  allows  one  to  recover  the  angular  dependence. 

We  employ  the  time  ordered  picture  of  ESD. in  implementing  this 
strategy.  An  ensemble  of  absorbed  ground  state  atoms  is  perturbed 
at  t  =  0  by  the  electron  beam.  The  forward  evolution  of  the  system 
is  followed,  without  further  interaction  with  the  external  beam,  to 
a  time  t  large  enough  for  all  ions  or  reneutralized  atoms  which  can 
desorb  to  have  crossed  a  detection  plane.  Each  trajectory  must  be 
individually  followed  to  its  asymptotic  direction  via  the  classical 
equations  of  motion.  One  evaluates  each  ion’s  or  atom's  contribution 
to  the  cross  section,  including  attenuation  due  to  quenching. 

Time  dependent  solutions  of  the  kinetic  equations  are  required. 
Hence,  we  obtain  formal  solutions  for  Pn  which  map  the  evolution 
of  an  initial  distribution  function  forward  in  time.  These  are  then 
used  to  construct  the  desorption  cross  sections. 

By  convention,  we  suppress  the  dependence  of  all  quantities  on  r, 
v,  displaying  only  their  parametric  time  dependence.  The  classical 
dynamics  problem  will  be  solved  numerically.  Liouville's  Theorem 


prompts  us  to  replace 

v-v+C(-~ 

D  -fc 


i, 
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after  which  they  take  the  form: 

57  = 

First,  we  introduce  the  ionic  survival  probability  against 
reneutralization,  which  is  defined  to  be: 
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The  formal  solution  pi(t)  for  the  ionic  state. is  found  to  be: 

r  * 

Kti)  =  p»  Xu)  s  jjco,  t)  +  R  ’  Ct.)  \  *t' c<’,< ) 

Jio 

Here,  we  integrated  the  equation  for  g(t)  formally  over  tQ^  t'  ^  t. 

The  boundary  condition  g(tQ)  =  P2(tQ)  ^oun<^  bY  taking 

the  limit  t  -»  tQ  +  £  .  One  notes  that  Lim  S2(tQ,t)  =  1.  Both 

t  ->  tQ 

R(t),  Q (t)  are  assumed  to  be  mathematically  well  behaved.  The 
equilibrium  distribution  satisfies  P°(t)  =  P°(tQ). 

Our  solution  for  P2(t)  yields  the  ionic  distribution 
function  which  has  evolved  from  the  initial  distributions 
P2(tQ)  and  p°(t0)*  The  first  term  displays  attenuation 
of  the  initial  ion  density,  while  the  second  incorporates  ionizations 
which  occur  during  the  interval  [vO- 

The  formal  solution  containing  the  ground  state  component  of 
quenched  particles,  is:  , 


f,  \t)  =  p/ft.)  + 1 P/CiOQ  Of1) 


Here,  one  uses  the  state  1  dynamics  to  map  r(t')r  v(t’)  forward  to 
time  t  in  the  second  term. 

We  commence  construction  of  the  cross  sections  by  first 


examining  the  response  of  a  single  ground  state  adatom  to  ionizing 
transitions  during  the  short  finite  interval  O-  °1  .  Let  U+ 

(— P , o)  be  the  probability  of  finding  the  atom  as  an  ion  at  t  =  0. 
This  is  found  to  be:  g 

U4('v)  ~  9-(v)  ^  ttr 0 ) 
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Here  we  chose  P^l-T)  =  0,  inasmuch  as  the  initial  ionic 
population  must  be  excluded  from  the  response  during  [-^3  .  w< 
normalized  to  unit  source  strength  with  the  choice  P°(tQ)  =  1. 

The  explicit  of  U  is  fixed  by  the  choice  of  initial  conditions 
r  (-*£}  ,  v(-T  ) . 

Next,  we  examine  the  evolution  of  a  single  ion  from  t  =  0  to 
time  t,  when  its  velocity  has  reached  the  asymptotic  value.  No 
further  interaction  with  the  external  electrons  is  permitted  to 
occur;  hence,  we  choose  P°(tQ  =  o)  =  o,  P^^o  =  °)  =  1* 

Conditional  probabilities  for  finding  the  atom  with  r(t),  v(t)  are 


given  by; 


V  5  C0/*) 


for  the  io  .ic  state,  and  £ 

V°M  - 

Jo 

for  the  ground  state.  An  ion  may  reneutralize  at  any  time  after 
t  =  o.  Hence,  V°(o,t)  contains  an  integration  over  intermediate 
times. 

o 

Next,  we  identify  single  atom  ejection  probabilities  T  ' 

(-t ,0,t)  in  the  ionic  and  neutral  state.  These  follow  the  history 
of  an  adsorbed  atom  which  responds  to  the  external  beam  over  M 

T  +(-T-  ,»,-{)=  <J*  (-■*,»)]/ +(o,<) 

-  {}+C-'tfo)v0  (0,i) 


These  quantities  are  closely  related  to  the  ionic  and  neutral  atom 

desorption  cross  sections  to  be  presented  below. 

New  we  introduce  the  ionization  cross  section  (t)  for  an 

o 

adsorbed  atom  at  the  point  r  =  r(t=0).  This  function  must  have  the 
free  atomic  or  molecular  ionization  cross  section  as  its  limit  when 
zq*^^  .  Here,  6  is  the  incident  electron  energy.  We  expect 

(rQ)  to  be  linear  in  the  upward  transition  rate  R  ^  as  follows: 

*«(£)  =  &)/ p 

This  result  is  recovered  after  first  noting  that  U+(-t,o)/t  is 
the  average  ionization  rate  per  atom.  The  corresponding  average 
ionization  cross  section  per  atom  is  given  by: 

*  J 

j. 

Now  the  limit  Tr*  o  is  taken,  to  obtain  the  instantaneous 
response.  R(t')  is  smoothly  varying  for  ESD,  inasmuch  as  the 
summations  over  final  state  electron  coordinates  contained  in  R^ 
consume  any  delta  function  resonances.  Ongfactors  R^(t'  =  0)  from 
the  integration  above,  and  makes  the  replacement  S2  (t '  ,0)  S2 
(-  X  ,0)  =  1. 

The  cross  section  for  the  entire  process  resulting  in  ionic 
desorption  of  an  atom  initially  in  the  ground  state  is  given  by: 

<rE;  M  =  CH.)  V+M  .  ■- 

Here,  we  explicitly  display  the  dependence  on  R  and  .Q.  The  step 
function  was  inserted  to  ensure  that  desorption  is  energetically 
possible.  One  integrates  the  decay  probability  appearing  in  the 
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exponent  from  point  to  point  as  the  trajectory  is  followed.  The 
asymptotic  velocity  generated  by  depends  on  the  choice  of  r (0) 

v(o).  One  chooses  t  for  each  ion  path  so  that  z(t)  lies  far  from 
the  surface. 

The  analogous  ground  state  ejection  cross  section  is  given  by: 

.  q  [*£*>  g*  fa)  -  st  (M + (<»n 

Here,  we  inserted  a  step  function  at  the  reneutralization  point  to 
eliminate  trapped  particles.  E2  governs  the  path  integrations 
during  (o,t')  ,  while  determines  the  dynamics  during  (t',tj  . 

Numerical  evaluation  of  (Tg7  is  expected  to  be  time  consuming  due 
to  the  integration  on  intermediate  times. 

Finally,  we  ensemble  average  the  ejection  cross  sections  over 
the  phase  space  coordinates  at  t  =  0.  The  formal  solutions  which 
yield  the  total  ionic  and  neutral  atom  ESD  cross  sections  are: 

<r + (e)  =  (a.)  ?,<v)  w] 

VCe)  P,°(a.^)7C 6^)  jfV  '  • 

v  0  0C/)  *  Still')-  Si  M"] 


The  equilibrium  distribution  P°  was  used  in  averaging,  rather 

than  P. ,  inasmuch  as  the  ejection  cross  sections  are  already 
1  + 
linear  in  R.  We  apply  our  expression  for  (T  below,  to  numerical 

simulation  of  ESDIAD  for  0+  adsorbed  on  Tungsten. 

Note  that  the  equation  for  ft'  is  a  generalization  of 

phenomenological  expressions  which  have  appeared  in  the  literature.. 
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Our  previous  work  on  one  dimensional  photodesorption  theory  produced 
analogous  forms  for  the  optically  fluorescent  and  ground  state  cross 
sections . 

The  result  expressed  above  is  model  independent,  except  for 
general  restrictions  on  R,  Q,  En,  and  P°  outlined  previously. 

One  may  use  the  above  equation  in  conjunction  with  realistic  model 
surface  interactions  to  obtain  ESDIAD  patterns  and  energy  spectra 
within  the  Menzel-Gomer ,  Redhead  ESD  picture. 

Finally,  we  note  that  the  transport  approach  is  applicable  to  a 
wide  variety  of  surface  phenomena  involving  propagation  of  atoms, 
ions  or  molecules.  In  particular,  the  cross  section  contributions 
to  higher  order  excitation  processes  in  ESD  may  be  obtained  without 
difficulty . 

Model  for  ionic  ESD  of  Oxygen  from  Tungsten 

We  solve  the  kinetic  equations  using  the  following  assumptions, 
for  simplicity.  These  do  not  change  the  qualitative  nature  of  our 
solution.  We  neglect,  for  the  time  being  all  quenching  effects,  and 
take  R  to  be  constant.  The  energy  curve  E^  is  chosen  to  be 

£<(£)=  -V0  4  £  (ft -l )  -K  •  £ ) 

where  are  fitted  to  best  experimental  "spring  constants".  The 

excited  state  energy  curve  E2(r)  has  an  attractive  (toward  the 
substrate)  image  potential  and  a  Hartree  potential  due  to  the 
tungsten  atoms  located  on  the  lattice  sites.  Thus  (in  atomic  units) 

) 


H  ^  K.  (&  ~  ~ 

A  ' 


dtK) 
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An  alternative  (and  opposite)  approach  to  the  ion  potential  is  given 
by  treating  the  tungsten  atoms,  or  lattice  positions  as  a  core  ions 
with  Z  =  6  and  letting  the  6  (2S,  4d)  electrons  provide  the 
screening  using  Thermes-Fermi  linearised  theory.  The  strength  and 
range  of  the  ionic  potential  are  about  the  same  as  obtained  using 
for  che  localized  tungsten  atomic  wave  functions. 

The  conclusion  of  our  calculations  indicates  a  large  focusing 
effect  of  ions  due  to  the  substrate  potential. 

The  first  phase  of  our  calculation  was  carried  out  for  the 
desorption  from  02/W.  In  our  phVs*  Rev*  Lett.  3j>,  610  (1976'  ,  we 
present  the  angular  distribution  curves  of  0+  ,from  W(100) 
surfaces.  We  thereafter  improve  our  computational  capability,  verified 
our  W(100) ,  results  and  obtained  new  results  for  the  angular  distri¬ 
bution  curves  of  0^  from  W(lll)  which  we  have  presented  in  Surface 
Science  6_9,  253  (1977).  The  most  attractive  feature  of  our 
calculations  has  been  to  show  from  first  principles  the  large  focusing 
effect  on  the  ions  which  results  in  sharp  angular  patterns  that  can 
be  traced  to  the  ion-solid  potential.  We  have  also  been  able  to 
show  that  for  geometries  where  strong  focusing  occurs  the  dominant 
part  of  the  repulsive  potential  is  due  to  the  nearest  metal  atom  or 
atoms.  One  difficulty,  however,  has  been  in  the  choice  of  the  inter¬ 
action  potential.  Due  to  the  lack  of  realistic  potential  curves  we 
have  adopted  the  simplest  possible  potential  imaginable  -  that 
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is  a  linear  combination  of  the  ground  state  atomic  Hartree  potential 
augumented  by  an  image  potential.  There  are  indications  that  this 
potential  may  not  differ  much  from  the  potential  surroundi-  /  a 
realistic  solid  and  thus  it  would  reasonably  represent  the  ion-solid 

4. 

interaction.  It  was  concluded  that  while  theory  predicts  for  0  from 

02/W  that  ions  should  emerge  predominantly  with  energies  of  2-4  eV, 

experiments  have  indicated  that  most  ions  are  ejected  typically  with 

energies  of  the  order  of  8  eV.  We  have  made  an  extensive  survey  of 

the  experimental  literature  and  have  found  that  the  maximum  kinetic 

energy  of  desorbed  0  from  02/Ni  an<^  °2/'Ge  about 

the  same  as  for  02/W  which  is  approximately  12  eV.  (See  our  paper 

Phys.  Rev.  16B ,  945  (1977)).  Repulsive  energy  of  12  eV  can  be 

++ 

accounted  for  by  the  existence  of  02  .  Although  our 

speculation  that  the  observed  0+  results  from  a  small  number  of 
non-dissociative  02  may  not  be  correct,  our  main  observation  that 

4. 

the  large  kinetic  energy  of  the  0  ions  may  result  from  doubly  (or 
triply)  ionized  "molecular  complex"  on  the  surface  is  yet  an  open 
question,  and  as  we  believe  may  yet  turn  out  to  be  true.  The  model 
of  Auger  induced  desorption  for  ionic  crystals  proposed  by  Knotek 
and  Fiebelman  (Phys.  Rev.  Lett.  -40,  964  (1978))  may  prove  to  be 
correct  also  for  covalent  as  well  as  ionic  bonding.  It  is  clear 
that  the  two-or  three-hole  intermediate  state  which  results  from  the 
Auger  process  of  core  ionization  for  covalent  bonding  will  be  long- 
lived  only  if  many  body  effects  are  dominant.  To  date  the  later  is 
also  an  open  question. 

4. 

Prompted  in  large  measure  by  our  finding  tv.r  0  desorbed  from 
02/W,  we  have  studied  the  ESD  of  H  from  H2/W.  Here  the 


%’• 
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metallic  characteristics  of  the  tungsten  substrate  are  preserved  and 
electronic  excitations  are  likely  to  disappear  in  times  shorter  than 
the  desorbing  time.  Our  results  (Surface  Science,  75,  L766  (1978)) 
for  the  energy  curves  of  H  from  Ii^/W  support  the  picture  of  ESD 
proposed  by  Gomer,  Menzel  and  Readhead. 

Our  research  has  widened  to  optic  and  electronic  properties  of 
surfaces.  We  have  calculated  the  interaction  of  polar-optical 
phonons  with  a  two  dimensional  electron  gas  (Surface  Science,  84 , 

440  (1979))  and  have  used  this  result  to  obtain  the  high  frequency 
conductivity  of  a  two-dimensional  electron  gas  interacting  with 
optical  phonons  (Phys.  Rev.  20B,  1306  (1979)).  This  result  is 
useful  for  "Metal-Oxide-Semiconductors"  devices  for  polar 
semiconductors.  We  heave  also  investigated  from  first  principles 
the  high-frequency  conductivity  of  a  two-dimensional,  two-component 
electron  gas  (Phys.  Rev.  20B,  4189  (1979)).  Here  we  have  considered 
the  effects  of  electron-electron  correlations  on  the  conductivity. 
Our  result  is  applicable  again  for  "Metal-Oxide-Semiconductor" 
devices  for  multi-valley  semiconductors. 

At  present  we  are  investigating  the  following  work:  (a)  the  high 
temperature  effects  on  surface  absorption  in  metallic  reflectors  for 
high  power  laser  use,  (b)  the  effect  of  energy  loss  of  charged 
particles  due  to  coupling  to  surface  vibrations  at  low  incident 
energies,  and  (c)  pumping  of  resonant  vibration  mode  of  admolecules 
on  a  substrate  for  very  short  pulses. 
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FIG.  1.  Sketch  of  the  (100)  face  of  tungsten  with  sev¬ 
eral  reference  points  labeled.  W  denotes  the  tungsten 
atoms. 


domly  selected  according  to  a  Maxwell-Boltznunn 
distribution.  The  atoms  are  converted  to  ions  (by 
the  excitation  process)  and  follow  an  appropriate 
dynamic  trajectory.  The  asymptotic  direction  of 
the  ion  is  recorded  and  displayed  in  the  figures 
to  follow.  For  those  cases  where  sharp  cones  oc¬ 
cur,  the  effect  of  surviving  the  quenching  proc¬ 
ess  just  leads  to  a  constant-numerical-factor 
correction.  Since  we  will  not  discuss  absolute 
cross  sections  here  but  only  the  angular  distri¬ 
butions,  we  will  neglect  such  numerical  factors. 
For  those  distributions  which  are  broader  in  ex¬ 
tent  the  effect  of  quenching  will  be  to  reduce 
sharply  the  flux  of  wide  angle  ions.  This  effect 
has  been  discussed  in  the  literature3  and  has  to 
do  with  the  amount  of  time  spent  in  the  vicinity 
of  the  surface. 


In  Fig.  1  we  present  a  sketch  of  the  (100)  sur¬ 
face  of  tungsten.  Various  points  on  the  surface 
have  been  labeled  for  reference  purposes.  Let  us 
start  by  examining  the  angular  patterns  produced 
by  assuming  a  single  adatom  lies  on  the  substrate. 
In  Fig.  2(a)  we  present  the  spot  pattern  produced 
by  assuming  that  the  original  atom  was  centered 
at  point  /I  at  a  distance  of  1  Bohr  from  the  sur¬ 
face.  Figure  2(b)  shows  the  spot  pattern  when  the 
atom  is  located  at  position  D  a  distance  of  4  Bohr 
from  the  surface  (both  positions  corresponding 
to  realistic0  W-0  bond  lengths).  We  see  that  the 
spot  patterns  do  not  have  lobes.  Positions  A  and 
D  are  high-symmetry  points  where  the  potential 
surfaces  have  approximate  cylindrical  symmetry. 
Thus  there  is  no  reason  for  preferred  azimuthal 
emission  directions  to  occur.  A  similar  result 
has  been  found  for  the  atom  at  point  E,  a  dis¬ 
tance  of  3.0  Bohr  above  the  surface.  In  the  low- 
coverage  phase  (/32  phase)  it  has  been  conjectured 
that  the  atoms  go  into  position  D  and  no  sharp 
lobes  have  been  found  experimentally,  although 
there  appears  to  be  a  hint  of  a  hazy  cross  pattern. 

In  Fig.  3(a)  we  present  the  spot  pattern  asso¬ 
ciated  with  the  atom  lying  originally  at  point  B 
(coordinates  .v  =  0.7</0,  .y  =  0.7uo,  z  =  1.0e/n,  where 
«=5.97«„  is  the  tungsten  lattice  constant).  Here 
we  see  sharp  islands  appear.  The  angular  pat¬ 
terns  have  a  fourfold  rotational  symmetry  and 
a  fourfold  reflection  symmetry — resulting  in  a 
replication  of  the  pattern  of  the  first  octant  in  the 
other  seven  octants.  In  Fig.  3(b)  the  spot  pattern 
associated  with  point  C  [x=  1.4«0,  y=0,  z=1.0ao) 
is  presented.  The  pattern  also  displays  sharp 
cones  but  rotated  by  an  angle  of  tj/4.  In  both 


FIG.  2.  Spot  patterns  corresponding  to  reference  points  A  and  0  of  Fig.  1. 
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FIG.  3.  Spot  patterns  corresponding  to  reference  points  U  and  C  of  Fig.  1. 


cases  the  upper-state  potential  surface  is  highly 
anisotropic  and  thus  the  dynamical  evolution  of 
the  trajectories  could  be  expected  to  produce 
sharp  angular  patterns.  The  degree  of  angular 
divergence  of  the  cones  as  well  as  their  size  de¬ 
pend  on  both  the  distance  away  from  the  surface 
and  the  distance  away  from  symmetry  points  in 
the  plane  of  the  surface.  Since  we  are  using  a 
rather  crude  ionic  potential  we  shall  not  explore 
the  systematics  of  the  variation  here.  However 
if  a  more  accurate  ion-solid  potential  were  avail¬ 
able,  one  should  be  able  to  pin  down  the  approxi¬ 
mate  location  of  the  adatom. 

One  must  emphasize  that  the  spot  patterns  are 
not  unique  to  the  geometries  cited  above.  If  one 
now  looks  at  situations  where  there  are  two  ad¬ 
atoms  (c.g.,  an  Oj  molecule)  adsorbed  on  the  sur¬ 
face,  it  is  possible  that  the  second  oxygen  atom 
provides  an  asymmetric  potential  for  the  io;.  to 


a 

FIG.  -t.  Spot  patterns  corresponding 


be  accelerated  by.  In  Fig.  4(a)  we  show  the  spot 
pattern  resulting  from  a  molecule  situated  at 
point  1)  with  orientation  parallel  to  the  [10]  direc¬ 
tion.  In  Fig.  4(b)  we  have  the  molecule  rotated 
by  an  angle  of  n/4  in  the  surface  plane.  The  inter¬ 
atomic  oxygen  distance  has  been  taken  to  be  that 
appropriate  to  molecular  oxygen.  Since  experi¬ 
mentally  the  cones  are  associated  with  the  f)t 
phase  of  adsorbed  oxygen  and  since  this  phase  is 
believed  to  involve  molecular  oxygen,  we  see  that 
a  reasonable  explanation  can  be  had  in  the  context 
of  such  a  model. 

In  summary  we  have  shown  that  the  angular  pat¬ 
terns  can  be  explained  in  terms  of  the  dynamics 
of  the  ion-solid  interaction.  Two  possibilities 
exist  to  explain  the  patterns.  Both  involve  the 
three-dimensional  spatial  asymmetry  of  the  po¬ 
tential  seen  by  the  oxygen  ion  when  it  is  produced. 
The  first  situation  involves  a  single  atom  located 
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at  a  nonsymmetry  site.  Such  a  nonconventional 
location  could  be  brought  about  by  subtle  long- 
range  interactions  between  oxygen  atoms  on  the 
surface.  The  second  situation  involves  an  oxygen 
molecule  pinned  to  a  tungsten  atom  on  the  surface 
The  present  theory  is  not  able  to  differentiate  be¬ 
tween  these  cases.  Thus  further  diagnostics  and 
analysis  are  necessary.  Currently  we  are  inves¬ 
tigating  the  possibility  that  vacancies  in  an  oxy¬ 
gen  overlaver  could  give  rise  to  additional  poten¬ 
tial  asymmetries. 
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A  first  principles  kinetic  theory  ts  developed,  applicable  to  angular  dependent  emission  in 
Electron  Stimulated  Desorption.  Expressions  for  the  ionic  and  neutral  atom  ESI)  cross  sections 
are  formulated  and  applied  in  a  model  calculation  of  O'  emission  from  W(l  1 1).  Strong  focus¬ 
sing  of  the  outgoing  ions  was  found,  with  the  use  of  a  model  ion-solid  potential  in  which  the 
substrate  was  free  of  excitation.  Off-axis  spot  groups  were  simulated.  The  peak  ion  energies  ob¬ 
tained  with  this  model  are,  however,  small  compared  to  experimental  energies  for  the  high 
coverage  case.  The  need  to  introduce  substrate  excitations  to  describe  this  case  is  discussed. 


1 .  Introduction 

Much  work  lias  been  done  in  recent  years  on  the  phenomenon  of  Electron  Stim¬ 
ulated  Desorption  (USD).  In  these  experiments,  a  beam  of  low  energy  electrons,  in 
the  range  of  100-500  cV  impinges  on  an  adsorbate  covered  surface,  causing  the 
ejection  of  ions  or  neutral  particles.  Several  articles  reviewing  ESD  have  appeared 
[1-4]. 

Experimental  data  have  been  collected  about  the  ion  yields  and  cross  sections. 
Energy  distributions  of  emerging  ions  have  been  measured  [5].  In  the  last  several 
years,  the  angular  distribution  of  ejected  ions  has  been  investigated  [6-8],  Strong 
focussing  of  the  desorbing  ions  has  been  revealed.  Desorption  patterns  were  ob¬ 
tained  which  contain  sets  of  well  defined  cones  for  ion  emission,  both  normal  and 
non-normal  to  the  surface.  The  number  and  orientation  of  emission  cones  is  related 
to  the  substrate  geometry. 

ESD  Ion  Angular  Distributions  may  thus  provide  valuable  clues  to  the  adsorp¬ 
tion  geometry,  on  the  scale  of  the  unit  cell  size.  This  direct  information  is  com¬ 
plementary  to  that  obtained,  for  example,  by  LEED,  which  is  sensitive  to  mainly 
long  range  order  in  the  adsorbed  layer  and  substrate. 

*  Research  Spnsuted  in  part  by  the  US  Office  of  Naval  Research  under  ONR  contract.  No. 

N000 14-75-00949  and  in  part  by  the  City  University  Faculty  Research  Award  Program. 
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In  this  theoretical  work,  we  address  ourselves  primarily  to  the  angular  depen¬ 
dence  problem.  Previous  theoretical  treatments  of  USD  -have  been  strictly  one¬ 
dimensional,  and  phenomenological  in. flavor  [9,10].  This  work  presents  a  first 
principles  transport  formulation  of  the  desorbing  ion  dynamics.  The  ion-solid 
interactions  are  assumed  to  be  time  independent,  or  long  lived  relative  to  ion  flight 
times.  > 

Section  2  establishes  the  conceptual  framework  and  introduces  kinetic  equations 
which  can  be  solved  formally,  in  three  dimensions,  for  the  linear  response  to  exter¬ 
nal  electron  bombardment.  In  section  3,  expressions  for  the  ionic  and  neutral  atom 
USD  cross  sections  are  obtained  from  the  kinetic  theory,  in  a  form  suitable  to  effi¬ 
cient  numerical  integration.  Section  4  is  an  illustrative  model  calculation  of  ESD 
Ion  Angular  Distributions.  The  experimentally  interesting  ease  of  0+  desorption 
from  W( 111)  was  chosen,  with  the  surface  assumed  free  of  excitation  during  the  ion 
ejection.  Section  5  is  reserved  for  discussion. 

We  found  that  ion  interactions  with  a  ground  state  surface  can  replicate  the 
focussing  and  ion  distributions  in  experimental  ESD  patterns,  for  several  adsorptiort 
geometries.  However,  this  model  surface  potential  probably  cannot  explain  the 
large  0+  ion  energies  observed  at  high  coverage.  We  speculatively  suggest  that  long 
lived  excitations  on  the  surface  may  play  the  dominent  role  in  ESD  of  0+  at  high 
coverage. 


2.  Formulation  of  the  problem 

This  calculation  employs  a  multi-step  picture  of  the  desorption  process  which 
has  been  found  useful  in  previous  ESD  investigations  [1-4,9-11].  Adsorbed  par¬ 
ticles  arc  Franek-Condon  excited  to  an  anti-bonding  state  from  which  ionic 
desorption  may  occur.  Propagation  follows,  with  the  dynamics  determined  by  an 
effective  adparticlc-surfacc  interaction  potential.  Recapture  may  occur  via  bond 
healing  transitions.  Rcncutralizcd  particles  may  still  possess  sufficient  kinetic 
energy  to  escape  as  neutrals,  in  which  ease  ground  state  propagation  completes  the 
desorption  process. 

The  above  picture  is  useful  for  describing  the  response  of  adparticle  substrate 
systems  which  are  interacting  weakly  with  an  external  beam,  lienee,  we  linearize 
the  kinetic  equations  presented  below  in  the  incident  flux.  Adparticle  propagation 
is  described  classically.  The  excitation  and  decay  mechanisms  involve  electronic 
transitions,  and  thus  require  quantum  mechanical  descriptions. 

A  single  atom  or  molecule  interacting  with  a  surface  may  be  described  by  the 
following  Hamiltonian: 

//=-  — +//c('V'/)-  (1> 

2  Al 

Here,  //e  contains  all  the  potentials  in  the  system,  as  well  as  kinetic  energy  opera- 
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tors  for  electrons  which  participate  in  atom-surface  interactions.  The  kinetic 
energy  opciutor  shown  acts  on  adparticle ’mass  center  coordinates,  with  At  denoting 
the  adparticle  mass.  V.‘  .  re¬ 

in  the  Born-Oppenheimer  approximation,  the  propagation  phase  of  desorption 
decouples,  to  order  ntJAl,  from  the  electronic  states  of  the  system,  liloctronic 
responses  to  shifts  in  the  adparticle  coordinates  may  he  regarded  as  instantaneous, 
permitting  a  set  of  effective  interaction  potentials  to  be  defined,  Excitation  and 
rcncutrali/.ation  arc  likewise  regarded  as  sudden,  inelastic  cs'ents  in  which  the 
internal  adparticle  state,  and  thus  the  propagation  dynamics,  change  discon- 
tinuously.  We  assume  that  the  asymptotically  recognizable  adparticle  slates  par¬ 
tially  diagonalize  lfe  close  to  the  surface. 

For  the  present,  consideration  is  restricted  to  a  two  level  system,  with  //  =  1  for 
the  neutral  ground  state.  For  positive  ions,  n  =  2.  The  z -dependent  behavior  of  the 


Fig.  1.  The  one-dimensional  qualitative  behavior  of  effective  potentials  for  ionic  ESD,  as  func¬ 
tions  of  r  coordinate.  1'0  is  the  chemisorption  energy.  The  r  coordinate  of  the  adsorption  site  is 
h.  The  metal  +  adsorbate  system  in  its  ground  state,  designated  by  (M  +  A),  has  the  potential 
function  A'j(r).  An  anti-bonding  potential  l:2{r)  leading  to  ionic  USD  is  denoted  (M  +  A+).  The 
substrate  here  is  free  of  excitation.  The  dashed  line  M*  +  A+)  suggests  the  behavior  of  the  effec¬ 
tive  potential  for  an  ion  interacting  with  the  excited  substrate.  The  asymptotic  energy  differ¬ 
ence  U2  -  A,  is  K|,  the  atomic  ionization  potential,  if  both  incident  and  ejected  electrons  arc  ’< 
Auger  emitted  to  the  vacuum. 
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effective  potentials  /:'„(/•)  is  sketched  qualitatively  in  fig.  1.  The  minimum  of 7;  j(r) 
lies  at  the  chemisorption  site.  In  ord^.r  for  ionic  desorption  to  occur,  the  antibond¬ 
ing  curve  /i2(r)  must  have  the  generally  repulsive  character  shown.  One  expects  IC2 
to  become  attractive  at  long  range  due  to  image  charge  density  residing  on  the  sur¬ 
face.  The  difference  (/sj(r)  --  /;  is  regarded  as  the  internal  ionization  energy  of 
the  adsorbate.  At  large  adsorbate  to  surface  separations  it  is  identical  to  the  free 
atomic  value.  Variations  of  /:'„(/•)  in  the  plane  paiallcl  to  the  surface  are  crucial  for 
the  angular  dependent  desorption  problem,  and  arc  discussed  below. 

This  picture  of  desorption,  and  the  transport  theory  developed  below,  may  be 
readily  genetalized  to  include  multiple  excited  state  systems,  in  which  substrate 
excitations  and  neutral  excited  states  of  adpartielcs  must  be  considered.  In  addi¬ 
tion,  photodcsorption  may  be  treated  within  basically  the  same  formalism  [12|. 

,  Little  attention  in  the  literature  has  been  given  to  the  ionic  states  which  govern 
the  desorption  kinetics.  Most  work  has  focussed  on  the  ground  state  potential 
minima,  where  adsorption  occurs.  Studies  of  electron  dissociation  of  molecules 
furnish  a  guide  to  the  antibonding  states,  however,  the  model  adopted  for  /f2  in  this 
preliminary  study  is  likely  to  be  crude.  The  difficulties  involved  in  modelling  l:2  for 
positively  ionized  cores,  arc  less  formidable  than  for  neutrals  or  negative  ions,  for 
which  exchange  and  correlation  effects  may  be  important. 

The  functions  /:’„(/■)  arc  assumed  known  in  developing  a  kinetic  theory  below. 
Furthermore,  the  solutions  of  the  classical  equations  of  motion  for  ions  and  neutral 
particles  arc  assumed  to  be  known,  and  to  adequately  describe  propagation  of  the 
massive  adpat  tides.  The  classical  solutions  r„(l)  depend  parametrically  on  time  and 
on  the  initial  cooidinates  for  each  adparticlc  trajectory. 

This  classical  description  is  most  applicable  to  unbound,  and  hence  desorbing 
states.  Ions  in  state  2  quickly  acquire  substantial  kinetic  energy.  The  same  argu¬ 
ment  applies  to  atoms  desorbing  in  the  ground  state  after  suffering  rcncutralization. 
A  classical  description  is  also  acceptable  for  vibrational  states  of  massive  adatoms  at 
high  temperatures,  that  is,  when  luo0  <  kT.  Were,  luo0  may  be  taken  to  be  the  level 
spacing  for  a  harmonic  oscillator  approximation  to  /;’(  near  the  chemisorption  site. 

When  /rw0  <:  kT,  classical  kinetic  equations  still  provide  qualitatively  correct 
results  for  the  linear  response.  The  propagation  of  particles  only  in  unbound  trajee^ 
lories  need  be  followed.  The  principal  correction  introduced  by  quantizing  the 
vibrational  states  is  to  modify  the  shape  of  ion  energy  distributions. 

We  denote  by  R(f),  Q(r)  the  transition  rates  for  ionization  and  rcncutralization 
respectively.  These  functions,  assumed  known,  of  adatom-coordinate  r  arc  averaged 
over  all  initial  and  final  states  of  the  substrate  and  incident  electron  consistent  with 
a  change  of  adatom  state.  Inasmuch  as  the  dominant  impulse  to  adsorbed' or  desorb¬ 
ing  particles  is  provided  in  this  picture  by  the  gradients  of  the  small  recoil  suf¬ 
fered  during  ionization  or  non-radiativc  decay  will  be  neglected. 

Electronic  energies  generally  employed  in  ESD  experiments  lie  in  the  range 
100-500  cV.  So  long  as  incident  electrons  arc  of  high  but  non-relativistic  energy, 
the  excitation  phase  of  ESD  may  be  described  in  the  Born  Approximation.  For 
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simple  ionization  piuccsses,  the  cooulinates  of  two  final  state  electrons  are  averaged 
over  in  R(r).  These  exit  quickly  from  the  surface  region,  and  arc  either  absorbed  by 
the  solid  or  emitted  as  Auger  electrons  to  thc-vacuum.  Inasmuch  as  incident  elec¬ 
tron  energies  arc  large  compared  to  (/;'2  -  E i),  ionization  can  occur  from  within  the 
entire  Franck- Condon  region. 

Non-radiative  decay  channels  dominate  transitions  back  to  the  ground  or  low 
lying  excited  states.  Optical  decays  arc  discounted,  as  their  lifetimes  are  long  com- 
pated  to  typical  ionic  flight  times,  ton  ^neutralization  involves  matrix  elements  for 
electron  capture  from  the  solid.  Hence,  this  process  is  efficient  at  close  range,  and 
weakens  rapidly  with  increasing  distance  from  the  surface. 

The  surface'  plane  is  assumed  to  be  tessellated  by  a  periodic  grid  of  of  surface 
unit  cells,  each  of  which  is  anisotropic  in  the  x—y  direction.  We  regard  this  struc¬ 
ture  as  known.  Facetting  of  the  suifacc,  or  rearrangements  of  the  substrate  layer 
spacing,  which  arc  known  to  occur  under  certain  conditions,  can  be  incorporated 
into  later  calculations  as  required. 

The  symmetry  and  periodicity  of  the  surface  mesh  arc  imposed  on  the  functions 
L'„,  R,  Q  which  describe  the  stages  of  desorption.  Transverse  anisotropies  in  the 
adparticlc-surfacc  interactions  become  ncgliblc  while  the  atom  to  surface  separa¬ 
tion  is  small  computed  to  the  sample  dimensions. 

We  formulaic  a  classical  transport  theory  entirely  within  a  single  surface  unit 
cell,  extended  in  the  z  direction  normal  to  the  surface.  At  macroscopic' distances, 
each  such  cell  behaves  as  a  point  source  of  angular  dependent  intensity.  Transport 
equations  linear  in  the  excitation  rate  R{r)  arc  now  obtained.  The  rate  Q(r)  is 
treated  as  a  zero  order  quantity.  Later  calculations  verify  that  the  desorption  cross 
sections  depend  exponentially  on  Q(r). 

Let  Pi,(r.v,t)  d3r  d3n  be  the  probability  of  finding  an  atom  or  ion  in  state  n  at 
time  /,  within  the  phase  space  volume  clement  at  r,  v ,  where  these  arc  the  ionic 
position  and  velocity.  The  following  normalization  is  assumed,  for  the  present: 

Tj  f  d3rd3u/,„(r,u,/)=  1  . 

"a'*2  unit 
cell 

This  condition  is  tclaxcd  after  linearizing  the  theory. 

All  spontaneous  desorption  processes  arc  neglected  for  the  systems  we  consider 
at  the  experimentally  interesting  temperatures.  Thermal  desorption  may  occur  as  a 
by-product  of  local  surface  heating  by  the  incident  electron  beam.  This  effect  may 
be  disicgarded  for  weak  fluxes. 

In  the  absence  of  external  perturbation,  the  atom-surface  system  is  taken  to  be 
a  gas  of  non-interacting  particles  in  thermal  equilibrium  with  the  solid.  State  2  is 
unpopulated,  hence: 

/1(r,  u)«0. 

Henceforth,  zero  superscripts  denote  equilibrium  state  quantities.  The- Boltzmann 
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equation  which  describes  the  ground  state  is: 


,/d/v 

Id/. 


4 


•0. 


(2) 


’  /SUIIK 

Here,  M  is  tlie  adsorbate  mass.  The  surface  interaction  term  describes  collisions 
between  adsorbate  atoms  and  substrate  excitations.  In  live  present  case,  these 
provide  principally  a  relaxation  mechanism  for  the  ground  state. 

Now  we  allow  the  system  to  be  driven  off  equilibrium  by  an  external  electron 
beam.  Competition  between  the  rclevcnt  transition  mechanisms  is  summarized  in 
the  following  kinetic  equations: 


^+u.Vo  ly,-. 
Dr  A~M  rL' 


-s_-a-a 


Here,  the  subscripts  Q.  E  refer  to  the  quenching,  excitation  mechanisms.  We  made 
the  choice  (d/^/dOsuiU'^O.  in  writing  the  above,  thus  neglecting  scattering  and 
energy  loss  of  the  desorbing  ions  due  to  interaction  with  substrate  excitations.  Such 
processes  will  be  studied  in  future  investigations. 

Ions  entering  state  2  are  described  by  the  source  term  (d/yd;)E.  Inasmuch  as 
recoil  is  neglected,  one  may  approximate: 


Tlie  corresponding  ground  state  depiction  occurring  via  the  ionization  process 
is  given  by: 

(d/»,/d0i:  =  i)  K(r)  ~ -(d/ydr),.; . 

Depletion  of  the  ionic  state  through  rcneutralizatic.i  is  given  by: 


(d/,:/d/)y  =  -/>2(M>,0  Q(r) 

Rencutralized  adparticles  which  re-enter  the  ground  state  arc  the  only  ones 
which  can  desorb  in  state  l.  These  contribute  a  distinctly  lion-thermal  component 
to  I ’|,  which  must  be  proportional  to  the  following  repopulation  term,  neglecting 
the  recoil  upon  quenching: 


(d/Vd/)y  *  P2{r,  u,  /)  (?(r)«  -(dF3/d/)Q 

The  following  set  of  kinetic  equations,  linear  in  R,  describes  the  angular  depen¬ 
dent  ESD  problem  for  the  weakly  perturbed  ease: 

D/V 


D  / 


+  «-•  W  -  ~  Vf/-2  •  VvP2’  =  P,°R  -  P2Q , 

/VI 


(3a) 


R.  Janow,  jV.  Tzoar  /  Kinetic  theory  description  of  USD  •  259 


9  ft.' 


a  t 


r*  »•  v,P,'--7  v;,  •  \j>;  =  P±Q. 


M 


(3b) 


All  terms  here  ate  of  order  ft1,  as  required.  One-dimensional  analogs  of  cqs.  (3a,  3b) 
were  employed  in  previous  work  on  photodcsorption. 

In  obtaining  the  result  above,  we  recalled  that  ft?®5 0.  Hence,  to  lowest  order 
ft2  sz  p'j  a  ft.  The  ground  state  distribution  was  taken  to  be: 

where  the  depleted  equilibrium  component  ft”  is  defined  to  satisfy: 


W£  +  tt.VP»  i  Vfvvri>=H 

I?  ;UV"  Vl  Id  th 


-PXR. 


SURF 


The  remaining  contribution,  /*j,  is  then  the  solution  to  cq.  (3b).  Inasmuch  as  ft” 
satisfies  an  equation  independent  of  <2,  a  complete  description  of  ground  state 
desorption  is  yielded  by  /y, .  Eq.  (2)  was  used  to  approximate  ft”  ~  ft1?  to  dominent 
order  in  ft.  For  weak  coupling  to  the  electron  flux,  the  characteristic  survival  time 
in  the  ground  state  against  ionization  is  long  compared  to  the  time  for  substrate 
induced  relaxation.  The  following  condition  then  holds: 


(dft,/df)su.u-Mdft,/dOi.:. 

Depletion  of  ft?  may  then  be  neglected,  as  its  effect  on  P'1  and  P\  is  of  order  ft2  at 
least. 

In  experimental  situations,  one  ideally  observes  the  current  of  ions  or  neutral 
atoms  ejected  in  each  state.  Linearity  of  ESD  ion  currents  with  the  incident  elec¬ 
tron  flux  has  been  supported  by  observation,  lienee,  the  intermediate  and  strong 
coupling  limits  of  this  theory  appear  to  be  of  little  practical  interest  at  present. 


3.  Desorption  cross  sections 

The  total  cross  section  for  desorption  in  the  nth  state  is  given  by: 
o„  =  /„(/..  t)IPQ(t) . 

Here,  /„(:./)  is  the  net  adparticlc  current  in  state  <  crossing  a  plane  parallel  to  the 
surface  at  z,  at  time  f,  per  unit  cell.  /„  is  evaluated  for  L  large  compared  to  unit  cell 
dimensions.  /„  is  proportional  to  both  the  electron  flux  /•’,  and  to  0(r),  the  number 
of  adsorbed  atoms  per  unit  cell  remaining  on  the  surface  at  t .  The  cross  sections  arc 
independent  of  /•’  for  linear  response.  The  cove-age  0(f)  and  current  I„(!-, l)  have 
the  same  time  dependence  via  hence,  o„  is  time  independent  as  well.  Here  a 
is  the  desorption  cross  section  summed  over  all  channels,  while  t  is  the  elapsed  time. 

For  angular  dependent  ESD  studies,  the  differential  desorption  currents /„(0,<ft> 
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r).  and  the  differential  cross  sections  o„(0,  0)  are  of  interest.  Mere  0,0  specify  the 
asymptotic  emission  directions.  Th'c  asymptotic  speeds  arc  integrated  over.  The 
total  and  differential  cross  sections -are  related  by; 

=f JJl  o„(0, 0) . 

One  may  find  explicit  expressions  for  o„(0, 0)  in  terms  of  the  distribution  func¬ 
tions  by  straightforwardly  following  the  microscopic  adparticlc  llux  within  one  unit 
cell: 


1 

o„(0, 0)  =  J*  f  duo3  f d \  Pn(rn ,7.  =  L u ) . 

‘  0  /1c 

Here  /■(  lies  m  the  plane  parallel  to  the  surface.  Ac  is  the  unit  cell  area.  One  first 
evaluates  the  current  density  contribution  from  d3u: 

J„(r,u,  r)  d3u  =  0(r)  /*,(r,u)  ud3u, 

associated  with  the  pitase  point  r,u. 

In  view  of  the  complexities  associated  with  the  angular  dependence  problem,  we 
adopt  an  approach  well  suited  to  numerical  simulation  of  USD  Ion  Angular  Distri¬ 
butions.  The  explicit  result  above  for  the  differential  cross  sections  is  cumbersome 
to  implement  efficiently  on  the  computer. 

The  steady  state  USD  response  is  found  by  solving  an  equivalent  time  dependent 
problem.  We  simulate  the  flux  of  ions  or  neutral  atoms  desorbed  in  response  to  an 
incident  electron  pulse  of  short  duration  r,  then  take  the  limit  r  •*  0.  Suitable 
expressions  for  the  total  desorption  cross  section*  are  constructed  using  the  trans¬ 
port  formulation.  Numerical  differentiation  in  the  asymptotic  plane  allows  one  to 
recover  the  angular  dependence. 

This  strategy  is  implemented  using  the  time  ordered  picture  of  USD.  An 
ensemble  of  adsorbed  ground  state  atoms  is  perturbed  at  t  =  0  by  the  electron 
beam.  The  forward  evolution  of  the  system  is  followed,  without  further  interaction 
with  the  external  beam,  to  a  time  r,  when  all  ions  or  reneutralized  atoms  which  can 
desorb  have  crossed  a  detection  plane.  Uach  trajectory  must  be  individually  fol¬ 
lowed  to  its  asymptotic  direction  via  the  classical  equations  of  motion.  One  eval¬ 
uates  the  contribution  made  by  each  ion  or  atom  to  the  cross  section,  including 
attenuation  due  to  quenching. 

Time  dependent  solutions  of  the  kinetic  equations  arc  required.  The  formal  solu¬ 
tions  for  P„  map  the  evolution  of  an  initial  distribution  function  forward  in  time. 
These  arc  used  to  construct  the  desorption  cross  sections. 

By  convention,  the  dependence  of  all  quantities  on  r,  o,  is  suppressed,  with  only 
their  parametric  time  dependence  displayed.  The  classical  dynamics  problem  will  be 
solved  numerically.  Liouvillcs  Theorem  prompts  the  replacement  in  eqs.  (3a),  (3b):  • 

3/dr  +  u  •  Vr  +  n  •  V„  -*■  d/dr . 
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Here. a  is  the  acceleration  in  the  nth, s,t'ate  given  by: 
I 

o  -  — n  V,  t:„(r). 

M 


Eqs.  (3a),  (3b)  now  lakr  the  form: 

tp'Am = i%)R(t)  -  UD  «o.  (4a) 

d/'I/tiz  =  />;(/)  Q(t).  .  (4b) 

First,  tlte  ionic  survival  probability  against  rencutralization,  S2(t o.  0*  >s  defined 
to  be: 


.Vj(/o.  0  ^  exp(-  J  dt'  GU')  1  •  ^ 

*0 

The  exponent  is  the  rencutralization  rate  integrated  along  the  ionic  trajectory  from 
r(/0)  to  r(r),  which  is  governed  by  f:2(r).  An  alternative  notation  which  displays  the 
path  integration  explicitly  is: 


•V2(r0,r)  =  cxp  -J 
L  '0 


C >(r')u(r') 
MO  I2 


Hcrc,u(r’)  is  the  ionic  velocity  at  r  . 

It  is  useful  to  interpret  ,V:(/0,  t)  as  an  evolution  operator,  which  maps  ions  from 
their  initial  to  final  coordinates.  One  obtains  the  form  of  cq.  (5)  by  solving  cq.  (4a) 
with  the  ionization  term  neglected. 

HSD  cross  sections  are  experimentally  found  to  be  several  orders  of  magnitude 
smaller  than  atomic  or  molecular  ionization  cross  sections  would  indicate,  lienee, 
quenching  is  presumed  to  be  efficient  and  we  expect  S2Uo>  0  I  •  when  r(t 0)  is 
close  to  the  surface,  and  r(t)  lies  in  the  detection  plane.  It  is  also  clear  that 
S2  (to.  0  O  ('.  t')  =  S2(t0.  /'),  and  S:(t0.  0  =  S2'(t.  t0). 

The  following  trial  solution  for  Pf2(t)  is  adopted: 


PkOadOS&o,/). 

After  formally  differentiating  the  above  and  comparing  with  cq.  (4a)  the  unknown 
function  g(t)  is  found  to  satisfy: 

dg(t')/dr'  =  /<(t')I/l(t,)X,(ro,r'). 

The  formal  solution  P'2(t)  for  the  ionic  stale  is 


P2(t)  =  P'2(to)S2(t0,t)  +  l*l(t0)  J  <lt'R{t')S2(t',i). 

>o 


(6a) 


Here,  we  integrated  the  equation  for  g(t)  formally  over  t0<t'  <t.  The  boundary 
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condition  g(to)~P’iU o)  >s  found  by  taking  the  limit  f-*t0  +  c.  One  notes  that 

lim  S2(to,  t)-  \  . 
fio 

Both  R(t),  Q(t)  arc  assumed  to  be  mathematically  well  behaved.  The  equilibrium 
distribution  satisfied  /^(/)  =  /' “(/„). 

The  solution  ly2( t)  yields  the  ionic  distribution  function  which  has  evolved  from 
the  initial  distributions  /’2(f 0)  and  /^(/o).  The  first  term  in  eq.  (6a)  attenuates  the 
initial  ion  density,  while  the  second  incorporates  ionizations  which  occur  during  the 
interval  |/0,  / ] . 

The  formal  solution  of  eq.  (4a),  containing  the  ground  state  component  of 
quenched  particles,  is: 


m*r\«o)*f4i'W)Qor 


Here,  one  uses  the  state  1  dynamics  to  map  /•(/),  u(r’)  forward  to  .time  t  in  the 
second  term. 

We  construct  the  cross  sections  after  first  examining  the  response  of  a  single 
ground  state  adatom  to  ionizing  transitions  during  the  short  finite  interval  [-r,  0]. 
let  (/(-  r,  0)  be  the  probability  of  finding  the  atom  as  an  ion  at  /  =  0.  Using  eq. 
(6a)  this  is  found  to  be: 


U\~tA))=  f  iU'/<(t')S2(t\0). 


Here,  the  choice  /y2(-r)  =  0  was  made,  inasmuch  as  the  initial  ionic  population 
must  be  excluded  from  the  response  during  |~r.  01.  We  normalized  to  unit  source 
strength  with  the  choice  /*?(/ o)  =  I  •  The  explicit  value  of  if  is  fixed  by  the  choice 
of  initial  conditions  r(  r),  u(-r). 

Next,  a  single  ion  is  permitted  to  evolve  from  (  =  0  to  time  t,  when  its  velocity 
has  reached  the  asymptotic  value.  No  furl h  r.interaction  with  the  external  electrons 
is  included;  hence,  we  choose  /’®(t0  =  0)  =  0,  and  P'£t0  =  0)  =  1  in  eqs.  (6a),  (6b). 
Conditional  probabilities  for  finding  the  atom  with  r(r),v(t)  are  given  by: 

C‘(0,/)  =  .V2(0,/) 

for  the  ionic  stale,  and 

i 

y°(o,n=f  dr's2(o,/W) 

o 

for  the  ground  state.  An  ion  may  rcncutra  izc  at  any  time  after  /  =  0.  Hence,  K°(0, 
t)  contains  an  integration  over  intermediate  times. 
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The  single  atom  ejection  probabilities,  7*,0(-t,  0,  t),  for  the  ionic  and  neutral 
state,  follow  the  history  of  an  adiOVbed  atom  which  responds  to  the  external  beam 
over  (-r,  Oj: 

7,‘(-r.0./)»t/,(-r.0)r',(0.0.  T\~t,0,i)=  U\-t,Q)  (8) 

These  quantities  are  closely  related  to  the  ionic  and  neutral  atom  desorption  cross 
sections  to  be  presented  bclosv. 

Lei  Xt(ro)  denote  the  ionization  cross  section  for  an  adsorbed  atom  at  the  point 
r0  =  r(/  =  0).  This  function  must  have  the  free  atomic  or  molecular  ionization  cross 
section  as  its  limit  when  r0  -*•  “.  Here,  e  is  the  incident  election  energy.  We  expect 
Xc(ro)  1°  bc  linear  in  the  upward  transition  rate  Rc  as  follows: 

Xc('o)  =  Ht(r0)H'  (9) 

Here,  /'is  the  incident  electron  flux. 

This  result  is  recovered  after  first  noting  that  U\-t,  0 )/r  is  the  average  ioniza¬ 
tion  rate  per  atom.  The  corresponding  average  ionization  cross  section  per  atom  is 
given  by: 

1  0 

X<(-r,°)  =  ~  /  d/7z<(/')Sj(/\0). 

-T 

In  the  limit  r-»0+  one  obtains  the  instantaneous  response.  /\c(/')  is  smoothly 
varying  for  USD.  inasmuch  as  the  summations  over  final  stale  electron  coordinates 
contained  in  R  consume  any  delta  function  resonances.  One  factors  /\’,(/'  =  0)  from' 
the  integration  above,  and  makes  the  replacement  .S2( t‘,  0)«.>j (-»},  0)=  1  to  ob¬ 
tain  cq.  (9). 

The  cross  section  for  the  entire  process  resulting  in  ionic  desorption  of  an  atom 
initially  in  the  ground  state  is  given  by: 

O|:j(0,0  =  Xcl>0)^(0,  0 

+  E2(r0)  -  EjH)  cxp[  -  / d/W)l  . 

o 

Here,  the  dependence  on  R  and  Q  is  displayed.  The  step  function  was  inserted  to 
ensure  that  desorption  is  energetically  possible.  One  integrates  the  decay  probabil¬ 
ity  appearing  in  tiic  exponent  from  point  to  point  as  the  trajectory  is  followed.  The 
asymptotic  velocity  generated  by  (jjij  depends  on  the  choice  ofr0,  u().  One  chooses 
t  for  each  ion  path  so  that  r(/)  lies  far  from  the  surface. 

The  analogous  ground  state  ejection  cross  section  is  given  by: 

4j(0,  0  =  x<(/o)  F°(0,  t)  =  XcO-o) /  d'W)  cxp[-J  dr"<2(r")l 

o  o 

X  0  (-}  Mvl  +  E2(ra)  -  i:2{r)  +  F.x{r)  -  £',(«»)  |  . 
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Here,  the  step  function  was  evaluated  at  the  reneutralization  point,  eliminating 
trapped  particles.  /:‘2  governs  the  path  integrations  during  (0,  f'],  while  E\  deter¬ 
mines  titc  dynamics  during  [/’,  r|.  Numerical  evaluation  of  o|’j  is  expected  to  be 
time  consuming  due  to  the  integration  on  intermediate  times. 

Finally,  the  ensemble  averages  of  ejection  cross  sections  over  the  phase  space 
coordinates  is  taken  at  /  =  0.  The  formal  solutions  which  yield  the  total  ionic  and 
neutral  atom  USD  cross  sections  are: 

ot  =f< dV„  i\  1%  o,  v0)  Xc(ro)-V2(0.  t)  Q[[Mo20  +  /f2(r„)  -  /:,(«)  | ,  (10a) 

t 

(7?  =  /d-Vo  d-’vu  l%th  v„)  xc(fo)  /  t')Q(t') 

0 

X  0|U/u$  +  /f2(r0)  -  l-2(r')  +  /f,(r’)  -  /;,(«)) .  (10b) 

The  equilibrium  distribution  Z’?  was  used  in  averaging,  rather  than  inasmuch  as 
the  ejection  cross  sections  arc  already  linear  in  R.  F.q.  (10a)  for  a  is  employed 
below  in  the  numerical  simulation  of  USDIAD  for  0+  adsorbed  on  tungsten. 

Note  that  cq.  (10a)  for  a*  is  a  generalization  of  phenomenological  expressions 
which  have  appeared  in  the  literature.  Our  previous  work  on  one  dimensional 
photodcsorption  theory  produced  expressions  analogous  to  Gqtjs  (10a,  10b)  for  the 
optically  fluorescent  and  ground  state  cross  sections. 

The  result  expressed  in  cqs.  (10a),  (10b)  is  model  independent,  except  for  gen¬ 
eral  restrictions  on  R,  Q.  E,„  and  1^  outlined  above.  One  may  use  cqs.  (10a),  (10b) 
in  conjunction  with  model  surface  interactions  to  obtain  USDIAD  patterns  and 
energy  spectra  within  the  Mcnzci-Gomcr,  Redhead  USD  picture. 

Finally,  note  that  the  transport  approach  is  applicable  to  a  wide  variety  of  sur¬ 
face  phenomena  involving  propagation  of  atoms,  ions  or  molecules.  In  particular, 
the  cross  section  contributions  due  to  higher  order  excitation  processes  in  USD  may 
be  obtained  without  difficulty. 

4.  Model  calculation  for  0+  desorption  from  tungsten 

The  remainder  of  this  paper  applies  the  first  principles  kinetic  theory  in  a  model 
calculation  of  desorption  from  tungsten.  We  attempt' to  understand  the  interesting 
angular  and  energy  dependent  USD  results  found  recently  on  W(lll)  [6]  by 
numerically  simulating  the  USD  patterns  and  ion  energy  spectra  to  be  expected  for 
atoms  adsorbed  at  various  sites  oh  the  surface.  Little  is  known  of  the  nature  of  pos¬ 
sible  antibonding  states  for  ions;  therefore,  simple  models  are  adopted  for  ion-solid 
and  atom-solid  interactions.  Absolute  cross  sections  were  not  calculated  due  to 
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lack  of  infoimation  about  wave  functions  on  surfaces.  We  concentrate  on  the  spec¬ 
troscopy  of  USD;  that  is,  in  this  paper  we  Itave  studied  only  the  polar  and  azimuthal 
angular  dependence  and  the  cnergyspectra  of  emerging  ions. 

Experimentally,  a  rich  diversity  of  ESD1AD  patterns  were  found  on  \V(1 1 1)  and 
W(100).  For  high  coverage,  these  consist  of  groups  of  sharp  Of  emission  spots, 
whose  geometry  is  a  function  of  adsoiption  and  annealing  temperature.  0+  emission 
normal  to  the  surface  was  found,  as  well  as  off  normal  spots,  whose  number  and 
orientation  dearly  relied  the  underlying  substrate  geometry. 

The  antibonding  state  potential  funnel  is  assumed  to  be' the  mechanism  which 
produces  focussed  spot  patterns.  For  a  number  of  the  observed  ESDIADs,  one  can 
discount  the  possibility  that  off  axis  emission  is  due  to  thermally  activated  sub¬ 
strate  facetting  [13].  Here,  only  flat  substrates  arc  considered.  ■ 

The  peaks  in  ion  energy  spectia  were  found  to  occur  at  8  to  9  cV,  on  both 
W(100)  and  W( 111)  [6-8],  This  agrees  with  earlier  studies  [5|  on  polycrystalline 
tungsten  which  report  a  spectrum  peaking  at  8.8  cV. 

The  qualitative  features  of  our  result  arc  unaltered  by  the  following  choice  of 
models.  Special  variation  in  the  excitation  and  rcncutrali/.ation  lates  was  neglected; 
thus,  we  chose  \(  and  S2  to  be  constants.  When  vibrational  states  arc  highly  local¬ 
ized,  the  distortion  of  o+  introduced  by  this  approximation  is  not  critical.  Ion 
energy  spectra  will  be  skewed  toward  lower  energies,  as  the  preferential  reneutra- 
lization  of  slower  ions  is  neglected. 

Near  a  possible  adsoiption  site,  we  modelled  /ij  as  a  parabolic  potential: 

/!',(/■)  »=  -  F0  +  \ (r  »  It)  K(r  -  l>). 

Here,  It  is  the  adsorption  site  coordinate.  The  chemisorption  energy  V0  does  not 
appear  explicitly  in  these  calculations.  The  spring  constant  tensor  K  was  taken  to  be 
isotropic,  for  simplicity.  One  fits  A'0(j  to  the  best  experimental  values. 

We  chose  for  !* j*  gaussian  distributions  both  in  position  and  velocity  spaces.  This 
choice  is  appropriate  for  particles  in  an  harmonic  potential  well  both  for  high  tem¬ 
perature.  and  in  the  zero  temperature  case.  The  shape  of  ion  energy  spectra  results 
from  a  folding  of  l*\  against  /:2  in  these  calculations.  However,  peak  ion  energies  arc 
insensitive  to  the  width  of  I This  sixfold  integration  on  phase  coordinates  r0„v0  in 
eq.  (10a)  was. accomplished  by  the  Monte  Carlo  method. 

In  modelling  /:2,  ions  were  approximated  as  point  charges  which  respond  to  the 
local  electrostatic  surface  potential.  The  character  of  the  dominant  anti-bonding 
state  in  USD  is  not  known:  for  simplicity,  we  chose  to  examine  one  which  leaves 
the  substrate  free  of  excitation.  /f2  is  taken  to  be  the  sum  of  the  unexcited  sub¬ 
strate  I  larticc  potential,  plus  an  image  contribution,  produced  by  averaged  .electron 
density  fluctuations,  viz: 

IW^qVmW+qltV^r-Ri). 

Ki 

Here,  r/  is  the  ionic  charge.  Bonding  effects  in  the  substrate  were  neglected.  The  stir* 
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fac.c  llurtrce  poteniiai  whs  represented  by  the  sum  of  free  tungsten  atom  Hartrcc 
potentials  centered  on  the  lattice  sites.  Similar  superpositions  have  been  the  starting 
point  for  self-consistent  band  structure  calculations  in  bulk  tungsten  [14—16). 

The  liciman-Skillman  wave  functions  [17]  were  used  to  evaluate. I'A  numeri¬ 
cally,  which  was  then  fitted  to  the  form: 

>a(t)  7.72  c-  ,-29r. 

All  quantities  are  in  atomic  units.  As  a  check,  we  evaluated  also  the  screened 
tungsten  atomic  potential  in  the  Fcrmi-Thomas  approximation.  At  realistic  W-0+ 
separations,  the  resulting  potential  is  somewhat  smaller  than  the  above, 

An  ion-surface  plasmon  hamil Ionian  [18,19]  was  used  to  numerically  evaluate 
the  image  potential  K(M,  for  a  number  of  values  of  2  andrs,  the  plasma  parameter. 
The  result  for  singly  ionized  atoms  was  fitted,  to  the  following  form: 

rj  «  -l/[4z  +  1.16  +  1 .34  rs  +  0(rs2)] . 

Image  charge  density  resides,  in  this  model,  on  an  approximate  electronic  boundary- 
of  the  metal.  Including  plasmon  dispersion  shifts  the  divergence  in  the  classical 
result  at  the  dielectric  surface  to  the  interior  of  the  metal.  ■ 

Some  consequences  of  this  ESD  model  on  W(1 1 1)  arc  now  examined.  Extensive 
numerical  computation  using  100  K  for  the  desorption  temperature  was  done,  at  a 
wide  variety  of  possible  adsorption  sites. 

The  most  striking  result  follows:  the  naive  model  constructed  for  the  anti-bond¬ 
ing  potential  strongly  focusses  ions  created  in  the  neighborhood  of  most  plausible 
candidates  for  adsorption  sites.  The  most  compact  simulated  O1,  emission  spots  are 
associated  with  the  largest  mean  ion  energies.  Long  range  image  attraction  tends  to 
defocus  the  spots.  However,  this  mechanism  is  competitive  with  focussing  by  the 
repulsive  component  of  /f2  only  for  very  low  energy  ions. 

Each  adsorption  site  candidate  produces  a  single  emission. spot.  For  higher 
energy  ions,  the  local  value  of  — V/s" 2  at  the  point  of  ionization  provides  the  dom* 
incut  impulse,  and  sets  both  the  central  eolatitude  and  azimuthal  orientation 
relative  to  the  substrate  for  the  resulting  emission  spot,  in  many  cases,  the  dom- 
incut  contribution  to  /f2  comes'  from  a  single  substrate  atom.  Emission  spot  geo¬ 
metry  then  reflects  the  chemisorption  bond  direction. 

Wc  were  successful  in  generating  simulated  ion  distributions  quite  similar  to 
those  which  were  observed  on  \V(l  11).  Success,  however,  in  replicating  the  visual 
features  of  ion  distributions  docs  not  guarantee  that  this  model  for  the  antibonding 
state  is  applicable.  In  particular,  the  simulated  ion  energy  spectra  which  were  pro¬ 
duced  peak  well  below  the  experimental  ion  energies.  This  was  the  case  even  for 
adsorption  site  candidates  whose  ESDIAD  pattern  simulations  closely  resemble 
those  observed.  This  finding  is  discussed  further  in  the  concluding  section;  For  now, 
wc  romai  k  that  remodeling  of  the  anti-bonding  state  is  required. 

Despite  the  limitations  of  this  model,  it  is  useful  to  explore  the  most  plausible 
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picture  of  high  coverage  oxygen  chemisorption  it  can  generate,  guided,  by  experi¬ 
ment.  (i* 

Whenever  the  anti-bonding  state,  potential  lias  the  qualitative  behavior,  of  our 
model,  the  following  rules  associate  the  gross  features  of.  ESDIAD  spot  patterns 
with  adsorption  geometry.  Emission  spots  directed  off  the  surface  normal  are  evi¬ 
dence  of  adsorption  at  a  point  where  /f2  is  assymmctric,  that  is.  away  from  the  one 
fold  degenerate  unit  cell  sites.  Here,  we  define-//  fold  degeneracy  to  mean  that  a 
particular  point  is  found  in  it  equivalent  locations  within  a  single  unit  cell  whose 
environments  are  identical  except  for  combinations  of  rotation,  reflection,  and 
translation.  For  example,  in  fig.  2  all  the  points  labelled  <*,0,7, 5  arc  degeneiate.  At 


pig.  2.  The  W(1 1 1)  surface  unit  cell,  large  spheres  depict  tungsten  atoms.  The  A  atoms  lie  in 
the  surface  layer.  I!  and  C  atoms  reside  in  the  second  and  third  planes  in  the  solid;  The  A-A 
distance  is  8.44  a0.  Dashed  lines  through  I),  C  are  bilateral  symmetry  axes.  "Primary”  ESD 
emission  lobes  for  0+  on  W(1 !  1)  have  been  observed,  with  their  substrate  registry  as  indicated 
by  the  trio  of  heavy  circles. 
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the  high  symmetry  (onefold  degenerate)  sites,  the  anti-bonding  potential  possesses 
approximate  azimuthal  symmetry,  and  thus  leads  to  emission  peaked  normal  to  the 
surface. 

One  would  reconsider  this  interpretation  if  long  lived,  highly  localized  substrate 
excitations  could  be  shown  to  dominate  in  the  anti-bonding  potential.  Such  excita¬ 
tions  might  be  localized  on  neighboring  adsorbate  or  substrate  atoms,  and  could 
result  in  non-normal  emission  of  atoms  adsorbed  at  symmetric  points. 

Within  the  model  we  used,  symmetrical  groupings  of  off  axis  spots  arc  explained. 
Assymmetric  points  in  the  unit  cell  arc  several-fold  degenerate.  An  adsorbate 
domain  on  the  surface  in  which  such  sites  are  occupied  should  be  found,  with  equal 
occupation  probability,  in  all  degenerate  orientations.  The  corrcsponding.ESDIAD 
patterns  should  thus  contain  spot  superpositions  suggestive  of  the  substrate  unit  cell 
structure. 

We  illustrate  with  applications  of  the  model  to  \V(l  1 1).  Fig.  2  shows  the  W(1 1 1) 
surface  unit  cell.  Some  of  the  lengths  pertinent  to  the  W-0  system  in  are  tabulated 
in  table  1 .  Here  we  are  guided  by  experiment.  Under  conditions  for  which  facetting 
of  the  surface  is  believed  not  to  occur,  trios  of  0+  emission  spots  were  found,  hav-. 
ing  the  substrate  orientation  shown  in  fig.  2.  These  were  found  in  conjunction  with 
normal  emission,  or  with  a  second  spot  trio,  rotated  180°. 

In  this  simulation,  adsorption  at  the  highly  assymmetric  points  such  as  a  or  (Tin 


Table  1 

Some  characteristic  lengths  in  the  tungsten-oxygen  system;  unless  noted  otherwise,  lengths  ate. 
taken  from  either  Slater  [20)  or  Panting  [21 ) 


Empirical  atomic  radii  (A) 

Oxygen  1 20) 

[21j 

Tungsten  (20) 

[21) 

bulk  tungsten  (bee)  lengths  (A) 
l.attiee  constant 
W-W  intcractomic  spacing 

W(l  1 1)  surface  lengths  (A) 
l.attiee  constant 
intcrplauar  spacing 

W(100)  surface  lengths  (A) 

Lattice  constant 
Intcrplanar  spacing 

W  ~0  bond  lengths  (A) 

Sum  of  empirical  atomic  radii  [20! 

1211 

Experimental  W-0  bond  length 
onW(ltl)  5.22 J 


2.08  *  0.07 
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fig.  2  can  produce  only  patterns  of  six  spots,  rotated  from  the  observed  orientation. 
Adsorption  directly  ovci  one  of  the,  highly  symmetric  points  A,  ii,  or  C  produces 
only  emission  in  a  narrow  cone  normal  to  the  surface.  Adsorption  elsewhere  along 
the  unit  cell  mirror  planes  is  expected  to  generate  spot  tiios,  with  primary  oriented  , 
spots,  as  in  fig.  2,  corresponding  to  adsorption  between  points  A  and  C,  while  180° 
rotated  trios  correspond  to  occupied  sites  between  C  and  B,  or  B  and  A. 


i 


-J- 

f 

I 


c 

Pig.  3.  Simulated  USD  spot  patterns  for  tow  coverage  of  oxygen  on  W(ltt).  (a)  shows  the  ion 
distributions  for  adsorption  1.95  A,  over  point  C  on  fig.  2;  (b)  corresponds  to  adsorption 
almost  equidistant  from  tungsten  atoms  at  A  and  C:  while  (e)  pertains  to  adsorption  between 
sites  C  and  I).  This  series  provides  the  best  fit  to  experimental  patterns  of  Madcy-ct  al. 
Seale  marks  correspond  to  5°  increments  in  asymptotic  eolatitude  0.  The  projections  =  tan  0 
cos  0, )’  =  tan  0  sin  C  was  used. 


270 


R.  Janow,  A'.  Tzoar  /  Kinetic  theory  description  of  USD- 


We  located  those  sites  which  .yield  the  best  lit  with  this  anti-bonding  potential, 
to  the  experimental  liSDIAD  patterns  and  ion  energies.  A  reasonable  lower  limit  on 
the  W-0  bond  length  was  chosen  as  1.95  A,  the  sum  of  Slater’s  atomic  radii 
[20-22].  The  colatitude  of  simulated  primary  lobes  was  fitted  to  the  experimental 
[6]  value  of  34°. 

The  most  plausible  site  candidate  found  for  the  origin  of  ‘primary’  spot  trios 
(i.c.,  those  with  the  orientation  shown  in  fig.  2)  lies  almost  equidistant  from  thd  A 
and  C  layer  tungsten  atoms,  in  fig.  2,  along  the  unit  cell  mirror  plane.  A  site  almost 
equidistant  from  the  C  and  B  layer  tungsten  atoms  is  this  models  best  choice  for  the 
generator  of  180°  rotated  spot  trios. 

For  botli  of  these  sites,  the  mean  energy  of  emitted  ions  was  calculated  to  lie  in 
the  range  2  to  2.5  cV,  which  is  significantly  smaller  than  the  experimental  value  of 
8—9  eV.  For  the  normal  emission  spot  produced  by  adsorption  directly  over  point 
C,  this  model  predicts  ions  energies  peaking  at  1.3  cV.  Note  that  these  calculations 
predict  peak  ion  energies  an  order  of  magnitude  below  experiment,  for  adsorption 
sites  directly  above  the  surface  layer  tungsten  atoms  at  A.  Fig.  3  shows  the  spot 
patterns  for  sites  mentioned. 


5.  Discussion,  conclusions 

The  first  principles  kinetic  theory  presented  above  for  angular  dependence  in  ESD 
provides  a  vehicle  for  testing  model  ion-solid  interactions.  Calculations  performed  us¬ 
ing  a  simple  model  for  the  ground  state  surface  anti-bonding  potential  predict  strong 
focussing  of  ions,  like  that  observed  in  JiSDIAD  experiments  involving  0+  emission; 
Models  of  0+  chemisorption  on  W  were  inferred,  by  fitting  to  the  experimental  pat¬ 
terns.  Unfortunately,  poor  agreement  between  our  calculated  ion  energy  distribu¬ 
tions  and  the  experimental  [5]  ones  weakens  the  plausibility  of  such  chemisorption 
models,  for  high  oxygen  coverage  on  tungsten.  However,  we  expect  these  model 
calculations  to  apply  to  low  oxygen  coverage  on  W  for  which  02  is  believed  to  dis¬ 
sociate.  It  is  likely  at  low  coverage,  that  the  metallic, substrate  will  support  only 
excitations  whose  lifetimes  arc  much  shorter  than  ion  ejection  times. 

Complementary  model  calculations  for  O1'  desorption  from  W(100)  were  done. 
Some  of  the  results  have  been  reported  in  the  literature  [1 1  ].  Here,  we  remark  that 
the  ion  energy  distributions  obtained  in  model  computations  peak  at  about  2.5  cV, 
as  for  \V(  111). 

Substantially  the  same  picture  of  adsorption  results  from  a  range  of  models  for 
the  antibonding  potentials  which  arc  time  independent,  and  which  involve  basi¬ 
cally  scale  and  range  changes.  More  careful  modelling  of  the  potential  for.  ground 
state  surfaces  may  moderately  increase  the  calculated  ion  energies.  Fo^  example,  we 
replaced  the  tungsten  wave  function  by  a  set  which  includes  L-b  ccurpijDg  [23], 
then  repeated  calculations  for  some  selected  sites,  ion  energies  were  increased  by  a 
factor  of  roughly  1 .5  to  2,  with  negliblc  change  in  the  simulated  ESDfAD'pattcrns. 


R.  Janow,  N.  Tzoar  /  Kinetic  theory  description  of  ESD 


271 


It  appears  unlikely,  however,  that  a  model  anti-bonding  state  in  which  the  sub¬ 
strate  and  remaining  adsorbate  atoms  arc  free  of  excitation  can  adequately  explain 
the  high  ion  energies  observed  for  oxygen. 

It  is  interesting  to  note  that  ion  energy  spectra  peaking  at  1.7  cV.  were  reported 
for  ESD  of  H+  from  tungsten  [5],  This  is  consistent  with  the  energy  range  predicted 
by  our  model. 

We  speculatively  suggest  that  excited  antibonding  states  involving  the  surface 
and  adsoibatc  atoms  may  be  able  to  explain  the  large  0*  ejection  energies  found  at 
high  covciage.  Explicit  detailed  treatment  of  these  as  yet  unknown  states  requires 
cluster  calculations  beyond  the  scope  of  this  work. 

The  type  of  processes  we  suggest  are  analogous  to  those  involved  in  electron  dis¬ 
sociation  of  free  molecules.  A  collision  with  an  incident  electron  leaves  the  local 
surface  in  a  state  from  which  an  adsorbate  ion  dissociates.  The  local  region  is  left  as 
an  excited  fragment,  which  produces  the  antibonding  effective  potential-governing 
ion  propagation.  Excited  states  may  be  distributed  over  an  extended  cluster,  or 
localized  to  one  or  several  neighboring  atoms. 

Surface  excitations  can  produce  ion  energies  in  the  range  observed  for  0+  on  W, 
provided  their  lifetime  is  long  compared  to  times  spent  near  the  surface  by  ions. 
Picture  a  simplistic  process,  for  example,  which  leaves  both  0*  and  W+  on  the  sur¬ 
face.  The  coulomb  interaction,  at  1.9S  A  separation,  yields  7-8  eV  energies.  For 
O2  molecules,  energies  are  ~12  cV. 

In  this  picture,  the  lifetimes  of  surface  excitation  must  be  long  compared  to 
ionic  flight  times.  Otherwise,  the  strong  anti-bonding  potential  shuts  off  immediat¬ 
ely  after  its  creation,  leaving  the  ground  state  surface  potential  to  provide  the 
dominent  impulse  to  desorbing  ions,  as  in  the  model  calculations  presented  above. 
Long  lifetimes  for  surface  excitations  may  be  justified  if  one  postulates  a  transition 
in  surface  behavior  from  metallic  to  insulating,  as  coverage  approaches  a  full  mono- 
layer  for  oxygen.  Note  that  the  experimental. 0+  energies  pertain  to  the  high  cover¬ 
age  I?,  phase. 

Insight  into  the  physical  reasons  for.a  low  yield  02  and  3  high  yield  (3|  oxygen 
phase  ca-t  be  gained,  if  this  insulating  barrier  picture  is  a  realistic  one  for  high 
coverage.  In  addition  to  supporting  long  lived  excitations,  such  a  barrier  would 
inhibit  ion  rcncutralization.  Both  higher  ion  energies  and  smaller  hopping  matrix 
elements  would  appear  in  the  rcncutralization  rate,  compared  to  a  low  coverage  sur¬ 
face.  Dramatic  increases  in  the  ion  survival  probability  would  result,  due  to  the 
exponential  dependence  of  .*>2  on  (7(z). 

ESD  emission  of  ilh  from  W  can  probably  be  explained  using  the  ground  state 
surface  model,  with  no  need  to  postulate  the  formation  of  an  insulating  barrier  at 
high  coverage.  The  H+  yield  for  both  (I,  and  02  hydrogen  phases  is  small,  and  com¬ 
parable  to  that  of  oxygen  in  the  low  coverage  phase. 

We  can  devise  speculative  mechanisms  which  produce  geometrical  spot  patterns 
in  ESDIADs,  whether  it  is  the  ground  state  surface  potential  or  an  excited  anti¬ 
bonding  state  which  governs  ion  emission.  If  the  anti-binding  state  is  associated 
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with  highly  locali'/cd  excitations,  off  axis  spots  can  result  even  in  the  case  where 
adsorption  occurs  at  the  highly ! symmetric  points  of  the  unit  cell.  Note,  for 
example,  on  fig.  2,  that  an  adsorbate  ion  created  above  point  C,  together  with  exci¬ 
tation  of  an  adjacent  A  layer  tungsten  atom,  can  produce  downward  facing  spot, 
trios.  If  excited  states  are  multi-centered,  the  anti-bonding' potential  will  have  essen¬ 
tially  the  same  symmetry  as  the  giound  state  potential.  Off-axis  spot  groups  then 
correspond  to  adsorption  at  points  in  the  unit  cell  away  from  the  symmetry  sites. 

The  speculative  nature  of  some  of  our  remarks  indicates  the  need  for  accelerated 
theoretical  study  of  USD  mechanisms,  after  which  unambiguous  adsorption  site 
determinations  using  HSD1AD  may  be  made. 
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Relation  between  electron-stimulated  desorption  and  nondissociative  chemisorption* 
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A  relation  is  fouml  between  the  maximum  ion  .energy  associated  with  electron-stimulated  desorption  (ESD) 
ami  the  chemisorption  bond  length.  This  relation  suggests  that  the  large  kinetic  energy  of  the  desorbed  ion 
may  lie  understood  in  icrmx  of  ESD  of  undissociated  molecules. 


The  phenomenon  of  electron-stimulated  desorp¬ 
tion  (ESD)  of  tons  from  solid-state  surfaces  has 
been  a  topic  of  considerable  interest  within  the  last 
decade. 1,2  The  adatom-substrate  system  is  as¬ 
sumed  to  be  promoted  by  the  incident  electron  to  a 
repulsive  energy  curve  resulting  in  the  ejection  of 
an  ion  from  the  surface.  The  ions  are  emitted  with 
a  rather  broad  energy  distribution  extending  up  to 
some  maximum  energy  Eau.  In  this  paper  we  will 
show  how  important  information  concerning  the 
chemisorption  bond  may  be  obtained  from  a  know¬ 
ledge  of  Em„. 

We  shall  be  concerned  with  molecular  gases  ab¬ 
sorbed  on  various  substrates.  For  some  systems 
it  is  believed  that  the  molecule  chemisorbs  as  a 
unit,  while  for  other  systems  there  is  evidence 
that  dissociation  occurs,  and  that  the  chemisorbed 
species  is  atomic,  at  least  at  low  coverages.  Most 
of  the  experimental  information  concerning  ESD 
which  we  shall  discuss  comes  from  high-coverage 
conditions,  where  the  bonding  situation  is  not  that 
well  understood. 

In  order  to  interpret  l:mu,  one  must  have  an  in¬ 
sight  into  the  nature  of  the  repulsive  states  involv¬ 
ed  in  ESD.  Before  considering  the  case  of  a  mole¬ 
cule  absorbed  on  a  surface,  consider  the  simpler 
problem  of  dissociative  ionization  of  gaseous  Hs. 

In  the  electron-impact  event,  both  of  the  molecular 
electrons  may  be  regarded  as  being  excited  to 
higher  orbitals.  There  is  clearly  a  rich  abundance 
of  excited  states' available  culminating  in  those 
states  where  the  molecule  is  left  doubly  ionized. 

In  most  of  the  excited  stales  the  electronic  wave 
function  expands  to  a  size  where  it  is  no  longer 
effective  in  screening  the  internuclear  repulsion. 

If  s  denotes  the  bond  length,  one  expects  ICmu  to 
be  roughly  given  by 

~c  /S  ,  (1) 

where  Enll  represents  the  upper  limit  to  the  kinetic 
energy  available  to  the  dissociating  atoms. 

The  situation  for  more  complicated  diatomic 
molecules  cannot  be  radically  different.  While  it 
is  possible  to  excite  the  more  tightly  bound  orbitals, 


the  cross  section  for  sue'-  a  process  will  be  neglig¬ 
ibly  small  compared  with  that  for  disruption  of  the 
highest-lying  orbital  pair.  Thus  one  expects  Eq. 

(1)  to  remain  valid. 

In  the  case  of  a  free  molecule,  the  repulsive 
energy  is  partitioned  between  the  dissociating  par¬ 
ticles  in  a  maimer  consistent  with  energy  and  mo¬ 
mentum  conservation.  When  the  molecule  is  ab¬ 
sorbed  on  a  surface,  the  situation  is  more  compli¬ 
cated.  Let  us  invoke  a  model  where  t.ne  ion  de¬ 
sorbs,  whereas  its  partner  atom  (or  ion)  remains 
bound  to  the  surface.  It  is  convenient,  but  not 
necessary,  to  think  of  this  binding  as  being  har¬ 
monic  in  nature.  As  the  ion  desorbs,  the  bound 
atom-solid  system  may  be  left  in  any  of  the  ex¬ 
cited  vibrational  states.  It  is  important  to  realize, 
however,  that  there  is  a  finite  probability  of  being 
left  in  the  ground  state,  i.c.,  that  recoilless  r/e- 
sorl>liou  may  occur.  This  implies  that  the  maxi¬ 
mum  energy  available  to  the  desorbing  ion  is,  in 
fact,  still  given  by  Eq.  (1).  There  is  no  recoil  cor¬ 
rection  that  has  to  be  made  to  this  energy. 

There  is  some  indirect  evidence,  in  fact,  that 
this  is  the  case.  If  one  studies  the  energy  distri¬ 
bution  of  the  emitted  ions  as  a  function  of  the  pri¬ 
mary  electron  energy  (Lozier  plot)  near  the  thres¬ 
hold,  one  experimentally  finds  a  slope  of  unity.5 
The  interpretation  given 5  is  that  the  mass  of  the 
recoiling  complex  is  very  large,  but  it  is  clear  that 
it  could  also  be  interpreted  as  a  manifestation  of 
recoillcss  desorption/1 

In  Table  I  we  have  assembled  experimental  data 
for  various  molecules  absorbed  on  several  sub¬ 
strates.  The  molecular  bond  lengths  S,  were  ob¬ 
tained  from  gas-phase  data.*  We  have  also  includ¬ 
ed  an  estimate  of  the  atom-substrate  bond  length, 
S.,  by  averaging  the  homonuclea.  diatomic  bond 
lengths.®  If  the  molecule  remains  undissociated, 
then  e2/Sl  is  the  appropriate  value  for  Hml„  where¬ 
as  if  it  docs  dissociate,  it  is  likely  that  c'/S,  is 
the  correct  value  (assuming  that  the  substrate 
atom  can  become  singly  ionized). 

Comparison  of  the  theoretical  predictions  with 
the  experimental  values  for  Emlt  suggest  that  in 
most  cases  quoted  here  the  results  may  be  under- 
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TABLE  I.  Comparison  of  theory  ami  experiment. 


atecule 

Substrate 

Reference 

(eV) 

s. 

(r»0) 

St 

w 

cVs, 

(eV) 

0  Vs, 

(eV) 

Dissociated? 

O, 

W 

7 

13 

2.26 

3.9 

11.9 

7.0 

No 

8 

12.5 

o2 

Mo 

1 

11.3 

2.28 

3.72 

11.9 

7.3 

Both 

8 

7.5 

o» 

NI(IIO) 

8 

11 

2.28 

3.50 

11.9 

7.7 

No 

Oi 

Si(l  11) 

9 

11.3 

2.28 

3.3'. 

11.9 

S.t 

No 

02 

Ge(l  It) 

•J 

11. G 

2.28 

3.4G 

11.9 

7.9 

No 

CO 

W 

7 

12 

2.12 

3.9 

12.8 

7.0 

No 

CO 

Graphite 

10 

12.2 

2.13 

2.13 

12.8 

12.8 

No 

»! 

W 

7 

7 

1.40 

3.3 

19.4 

3.3 

Yes 

NO 

W 

11 

11.2 

2.18 

3.9 

12.5 

7.0 

No 

11,0 

W 

7 

15 

1.81 

3.3 

15.0 

8.3 

No 

CO, 

w 

7 

10 

2.53 

3.9 

10.8 

7.0 

No 

stood  in  terms  of  ESD  due  to  undissociatod  mole¬ 
cules.  Exceptions  to  this  occur  for  lla  on  tungsten, 
where  reasonable  agreement  with  experiment  is 
obtained  for  the  dissociation  model.  The  case 
of  02  on  Mo  is  interesting  because  it  appears  to  ex¬ 
hibit  both  dissociated  and  undissociatcd  molecules 
under  different  experimental  conditions. 

The  model  presented  here  is  extremely  simple. 

It  neglects  effects  due  to  image  forces  and  residual 
surface  fields.  The  fact  that  the  model  works  so 
well  can  be  taken  as  evidence  that  these  effects  are 
not  Important.  This  is  consistent  with  the  fact  that 


the  surface  is  strongly  oxidized  in  the  high-cover- 
age  pluses  of  chemisorption. 

In  the  present  discussion  we  are  neglecting  the 
smearing  of  the  high-energy  cutoff  due  to  the  zero- 
point  fluctuations  in  the  position  of  the  ground  state 
adatom.  This  wilt  be  small  due  to  the  slow  spatial 
variation  of  the  Coulomb  potential.  The  full  width 
at  half  maximum  for  this  is  estimated  to  be  AE 
=  (|iwsVlnlG)‘l/s,  where  p  is  the  reduced  mass  of 
the  adatom,  anti  w  is  the  vibrational  frequency  of 
the  undissociated  molecule.  For  CO  and  Oa,  for 
example,  A£  is  0.9  eV. 
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Electron  Stimulated  Desorption  is  a  phenomenon  in  which  a  beam  of  low  energy 
electrons  (100-500  cV)  impinges  on  an  adsorbate  covered  surface,  causing  the  ejec¬ 
tion  of  ions  or  neutral  particles.  Ion  yields,  ESD  total  cross  sections,  and  energy 
distributions  have  been  found  experimentally  for  a  variety  of  adsorbate-substrate 
combinations.  Most  recently,  angular  dependent  patterns  of  strongly  focussed  ion 
impact  spots  have  been  observed,  whose  geometry  may  reflect  the  substrate  unit 
cell  structure  [  1  — 3  J .  ESD  data  provide  important  clues  to  the  chemisorption  geom¬ 
etry  and  insight  into  the  character  of  surface  interactions.  Several  recent  review 
articles  exist  (4-7). 

A  phenomenological  model  of  the  ESD  process  in  one  dimension  was  suggested 
(8,9)  by  Comer,  Mcn/el  and  Redhead  (GMR).  This  multistcp  model  pictures  absor- 
bale  atoms  to  be  promoted  via  Franck-Condon  transitions  from  a  bound  vibra¬ 
tional  level  to  an  excited  anti-bonding  state,  which  may  dissociate  with  the  emission 
of  an  ion. 

This  paper  presents  a  first  principles  calculation  of  the  kinetic  energy  distribu¬ 
tion  for  H*  ions  desorbed  from  tungsten  surfaces  via  ESD.  Here  we  amplify  previ¬ 
ous  work,  which  focussed  primarily  on  the  angular  dependence  in  0+  emission  [10, 
1 1],  A  kinetic  theory  applying  to  ESD  in  three  dimensions  was  formulated  in  our 
cailier  work  from  first  principles  using  a  generalisation  of  the  GMR  picture.  Ions 
were  assumed  to  move  classically,  with  their  dynamics  governed  by  an  adparticle- 
surface  effective  potential  Bond  healing  .cneutralization  mechanisms  may  cause 
ions  to  be  recaptured.  These  may  desorb  as  neutrals,  if  their  history  as  ions  provides 
sufficient  kinetic  energy.  The  electronic  transitions  involved  in  ionization  and 
quenching  (rcneutrali/ation)  necessitate  quantum  mechanical  descriptions  of  these 
processes. 

This  study  of  1I2/W  is  piomptcd  by  our  recent  findings  for  0+  desorbed  from 
02/W.  We  calculated  ESD  Ion  Angular  Distributions  and  energy  spectra,  using  sim¬ 
ple  moJels  for  the  interaction  of  0+  with  the  surface.  Despite  the  crudity  of  the 
model  potential  employed,  rough  agreement  of  calculated  peak  ion  energies  with 
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experiment  should  have,  been  obtained,  if  a  ground  state  surface  generates  the 
potential  rclcvcnt  to  emission  dynamics. 

Experimental  determination  of  the  stale  (high  coverage)  0+  energy  peak 
yields  8.8  eV,  and  a  large  spectral  width  which  constrasts  sharply  with  the  predic¬ 
tion  of  our  model  for  0+.  Strong  focussing  effects  and  the  main  features  observed 
in  ESDIAD’s  were,  however,  replicated  in  our  earlier  work  for  several  speculative 
oxygen  adsorption  geometries.  We  arc  currently  investigating  the  hypothesis  that 
either  molecular  absorption  or  long  lived  excited  states  on  the  surface  may  explain 
the  BSD  spectrum  of  0+. 

Note  that  the  experimental  energy  distribution  for  ESD  of  H*  from  H2/W  peaks 
at  1.7  eV  for  high  coverage  on  polycrystallinc  material  [12].  This  lies  in  the  range 
expected  with  a  ground  state  surface  model.  Attention  is  therefore  restricted  here 
to  calculating  the  H+  spectrum,  inasmuch  as  it  provides  a  critical  test  of  the  GMR 
picture  of  ESD  employed  in  previous  work.  The  detailed  calculations  closely  paral¬ 
lel  those  discussed  in  ref.  [10].  Reasonable’ values  were  adopted  for  some  of  the 
phenomenological  parameters  required. 

A  three-dimensional  model  for  the  anti-bonding  H+-W  surface  potential  was 
constructed  for  W(100),  W( 111),  and  W(110),  in  the  manner  described  in  earlier 
work  [10].  Here  H+  ions  arc  assumed  to  respond  to  a  time  independent  potential 
field  which  we  approximate  by  the  sum  of  unperturbed  tungsten  atom  Hartree 
potentials,  plus  an  image  contribution.  The  corresponding  forte  governs  ionic  mo¬ 
tion  following  dissociation  from  the  surface. 

Experimental  e\idcncc  indicates  that  H2  adsorbs  on  W  in  atomic  form  -  with 
recent  experimental  confirmation.  At  high  coverage  (~2  atoms/unit  cell)  the  char¬ 
acteristic  vibration  energy  is  135  MeV,  contrasting  sharply  with  that  expected  for 
molecular  vibrations. 

We  approximated  the  ground  state  potential  for  hydrogen  atoms  by  three- 
dimensional  parabolic  wells  at  each  of  several  test  adsorption  sites,  using  spring 
constants  matched  to  the  experimentally  determined  vibration  frequencies. 

At  300  K,  the  ground  state  vibration  level  spacing  is  large  compared  to  the 
thermal  energy.  Hence,  absorbed  hydrogen  atoms  arc  assumed  to  vibrate  in  the  zero 
point  harmonic  oscillator  state.  H  atoms  ionized  by  electron  impact  arc  represented 
by  a  minimum  uncertainly  wave  packet  corresponding  to  a  classical  Gaussian 
distribution  of  the  form: 

3 

P?(r,u)=n  {(2np,-i'/)~1  exp(-u//2f/) exp  {--fa -  bi)2/2pf} )  .  (1) 

t=t 

Here  b  locates  the  adsorption  sites.  The  widths  in  atomic  units  are  given  by 

Pi  =  (2 ;»«/)“ 1/2 ,  vt  =  (co//2m) 1/2 .  (2) 

The  parameters  used  in  this  calculation  yield  p/ «  0.24  a0.  The  finite  width  in  velo¬ 
city  space  was  included  but  lias  negligible  effect. 
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Fig.  1.  The  qualitative  behavior  of  the  effective  potentials  for  ionic  USD,  as  functions  of  the; 
coordinate.  VQ  is  the  chemisorption  energy.  The  z  coordinate  of  the  adsorption  site  is  b.  The 
metal  +  adsorbate  system  in  its  ground  state  is  designated  by  (M  +  A).  An  anti-bonding  poten¬ 
tial  leading  to  ionic  USD  is  denoted  (M  +  A+),  The  substrate  here  is  free  of  excitation.  The 
dashed  line  (M*  +  A+)  suggests  the  behavior  of  the  effective  potential  for  an  ion  interacting 
with  the  excited  substrate.  The  asymptotic  energy  difference  is  the  atomic  ionization  poten¬ 
tial,  if  both  incident  and  ejected  electrons  arc  Auger  emitted  to  the  vacuum. 


We  limited  our  calculation  to  the  shape  of  ion  energy  distributions,  which  arc 
obtained  by  folding  the  ion-solid  potential  E2  against  the  ground  state  distribu¬ 
tion.  Fig.  1  provides  a  one-dimensional  illustration,  with  gaussian  shape  of  P°x 
skewed  to  form  a  high  energy  tail  in  the  spectrum,  due  to  variation  in  the  slope  of 
E-j, .  The  spectral  shape  will  also  be  modulated  by  the  ion  survival  probability  against 
rcneutralization.  In  the  present  calculation,  reneutralization  effects  arc  taken  to  be 
independent  of  initial  position,  as  is  the  ionization  cross  section.  The  parameters 
required  to  construct  realistic  models  of  rcneutralization  are  not  yet  known. 

This  (GMR)  interpretation  of  kinetic  energy  spectra  is  a  good  description  of  ESD 
only  if  the  substrate  remains  in  its  ground  state  during  ion  propagation,  which  is 
assumed  here,  or  if  the  potential  is  produced  by  substrate  excitations  whose  life¬ 
times  arc  long  compared  to  ion  flight  times.  If  the  cross  section  to  produce  surface 
excitations  is  large  enough,  with  the  decay  lifetime  comparable  to  the  time  of 
flight,  another  width  mechanism  may  contribute  via  competition  between  the 
relaxation  and  rcneutralization  lifetimes.  Inelastic  recoil  involving  substrate  and 
adsorbate  atoms  may  also  widen  the  ion  spectra.  This  mechanism,  however,  should 
not  be  significant  for  light  ions. 

It  was  necessary  to  calculate  H+  desorption  spectra  for  each  of  the  W(100), 
W(l  1 1),  and  W(1 10)  faces,  inasmuch  as  the  polycrystalline  sample  used  by  Nishi- 
jima  and  Propst  [12]  had  an  unknown  mixture  of  faces.  We  calculate  spectra  for 
hydrogen  atoms  adsorbed  at  each  of  several  highly  symmetric  points  on  each  unit 
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Pig.  2.  (a),  (b)  and  (c):  calculated  ion  energy  distributions  for  I5SD  of  H+  from  li*  on  W(100), 
W(1 1 1),  and  W(1 10).  All  spectra  arc  normalized  to  unity,  with  energies  given  in  electron  volts. 
The  broken  curve  shows  experimental  results  of  Nishijima  and  I’ropst  [12].  On  each  surface, 
highly  symmetric  sites  were  chosen  for  calculation,  as  shown  by  the  Insets  for  each  surface 
geometry,  with  3.3  a0  taken  to  be  the  W-II  bond  length.  On  W(100),  site  A'  is  five-fold  coordi¬ 
nated,  site  E  is  a  two-fold  coordinated  bridge  site,  and  sites  A,  D  arc  singly  coordinated.  On 
\V(111),  all  sites  are  singly  coordinated  3.3 eg  above  W  atoms. On  W(110),  sites  B,  Care  bridge 
sites,  site  0  is  three-fold  coordinated  to  surface  layer  W  atoms  and  site  A  is  singly  coordinated 
over  the  surface. 
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chemisorption  site  coordinate, 
n  a  normalization  constant,  the  resulting  spectrum  [12]  is 
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Fig.  3.  Calculated  spectra  for  USD  of  D+  from  tungsten.  Here,  spectra  for  the  sites  directly  over 
surface  tungsten  atoms  are  shown. 


tralized.  On  the  other  hand,  fast  ions  are  those  created  closest  to  the  surface,  and 
which  interact  most  strongly  with  substrate  electrons  increasing  the  quenching 
probability.  The  sense  in  which  the  resulting  energy  spectrum  is  skewed  is  strongly 
model  dependent.  For  R  =  1/2,  quenching  is  uniform  but  for  R  >  1/2  the  relative 
intensity  of  spectral  tails  present  in  fig.  2  sould  be  reduced.  The  opposite  effect 
occurs  for  long  range  quenching,  defined  here  by  validity  of  the  condition  R  > 
1/2. 

Experimental  ESD  spectra  for  D+  from  D2/\V  might  be  useful  in  establishing 
rcncutrali/ation  effects,  due  to  mass  dependence  in  the  exponent  above.  The  vibra¬ 
tional  states,  however,  possess  complicating  mass  dependence.  The  prominent  tails 
in  some  of  the  calculated  H*  spectra  in  fig.  3  can  be  attributed  to  the  large  ampli¬ 
tude  of  zero  point  motion  for  light  particles.  Fig.  3  shows  calculated  D+  spectra, 
which  show  considerable  nartowing  compared  to  their  counterparts  in  fig.  2. 

Note  that  higher  order  processes  which  involve  surface  excitations  are  possible 
contributors  to  the  H*  spectrum  in  the  region  near  its  cutoff  above  7  eV.  The  limit¬ 
ing  energy  is  given  by  the  H+-W+  Coulomb  repulsion,  which  is  about  8  cV  for  3.3 
a0  separation.  Somewhat  lower  energies  are  expected  for  high  rydberg  surface 
states.  Some  preliminary  calculations,  using  a  delta  function  model  for  the  vibra¬ 
tional  state  distribution,  show  that  finite  width  spectra  can  be  produced  via  com¬ 
petition  between  the  surface  state  lifetime  and  the  rcncutrali/ation  rate.  A  site  such 
as  A'  might  produce  a  broadened  spectrum  via  such  processes.  Another  process  pos¬ 
sibly  contributing  to  the  broadening  makes  use  of  coupling  between  translational 
states  of  the  ions  and  vibrations  in  the  surface.  Significant  energy  loss  would  result 
from  the  creation  of  such  multiphonon  final  states. 

To  summarize,  these  results  for  ESD  of  11*  from  H2/W  support  the  theoretical 
picture  of  ESD  proposed  by  GMR  and  formalized  in  our  previous  work.  Despite 
the  crudity  of  the  models  used,  this  calculation  produces  energy  spectra  consistent 
with  experiment  for  several  surfaces  and  plausible  atomic  adsorption  sites.  In  par- 
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ticular,  our  model  H+— W  potential,  in  which  the  surface  remains  free  of  excita¬ 
tion,  provides  a  qualitatively  correct  description  for  ESD  of  H+  from  H2/W.  The 
lack  of  experimental  data  about  ESD  cross  sections  and  rcncutralization  effects  pro¬ 
hibits  us  from  making  definitive  statements  which  would  select  the  dominant  sur¬ 
faces  and  adsorption  sites  contributing  to  the  experimental  H+  spectrum. 

Richard  JANOW  and  Narkis  TZOAR 
Department  of  Physics,  The  City  College  of  the  City  University  of  New  York,  New 
York,  New  York  10031,  USA 
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The  interaction  Hamiltonian  for  a  two-dimensional  electron  gas  with  optical  phonons  was 
obtained,  for  an  electron  gas  on  the  interface  of  a  polar  semiconductor  in  a  metal— insulator- 
semiconductor  device.  We  found  that  the  electron-phonon  interaction  is  proportional  to  the 
inverse  square  root  of  the  phonon  wave  number.  The  dependence  of  the  electron-phonon  cou¬ 
pling  constant  on  the  dielectric  constant  of  the  insulator  is  determined. 


1 .  Introduction 

When  an  electric  field  is  applied  to  the  surface  of  a  semiconductor,  as  occurs  in 
metal- insulator-semiconductor  devices  under  bias,  a  strong  accumulation  of  elec¬ 
trons  on  the  surface  of  the  semiconductor  takes  place,  known  as  inversion  layer. 
These  electrons  arc  trapped  by  a  surface  potential  and  are  bound  in  the  direction 
normal  to  the  surface,  but  arc  free  to  move  on  the  plane  defined  by  the  semicon¬ 
ductor-insulator  interface.  Since  the  average  distance  between  electrons  within  the 
plane  is  larger  than  their  spread  normal  to  the  plane,  one  may  describe  the  behavi¬ 
our  of  the  electrons  in  the  inversion  layer  as  two-dimensional. 

The  study  of  .transport  and  optical  properties  of  two-dimensional  electron  gas 
is  an  active  area  of  research  »n  previous  work  mostly  non-polar  semiconductors 
have  been  considered  [1].  Thus  electron-impurity,  clectron-accoustic  phonons 
and  multivallcys  electron-electron  collisions  have  been  studied  [2-8).  However, a 
recent  interest  in  polar  semiconductors  motivates  the  present  paper  [9).  Here  an 
important  scattering  mechanism  results  from  the  interaction  of  longitudinal  polar 
phonons  with  the  electrons.  We  present  a  dielectric  formulation  [10-12)  for  the 
interaction  of  a  two-dimensional  electron  gas  and  polar  phonons.  We  derive,  using 
general  principles,  both  the  matrix  element  for  the  electron-phonon  scattering  and 
the  surface  polar  phonon  frequencies. 


*  Research  sponsored  in  part  by  the  US  Office  of  Nav.il  Research  under  ONR  grant  No. 
N00014-75-C-0949  and  by  The  City  University  Faculty  Award  Program. 
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2.  Excitations  and  interactions  in  semi-infinite  dielectric 

The  long  wavelength  excitations  of  a  semi-infinite  polar  crystal  can  be  charac¬ 
terized  by  a  scalar  potential.  Here  we  neglect  entirely  the  small  region  of  very  long 
wavelengths,  where  phase  velocities  of  the  excitation  roughly  equal  the  speed  of 
light,  and  retardation  effects  become  important.  The  excitation  spectrum  of  semi- 
infinite  crystals  and  their  interactions  with  electrons  is  describable  in  terms  of  the 
bulk  dielectric  function  c(k,  to).  We  think  of  these  excitations  as  potential  waves 
propagating  in  the  medium  and  inducing  polarization  charges. 

We  start  with  the  well  known  case  of  the  infinite  crystal.  Here  the  excitation 
spectrum  is  given  in  terms  of  the  longitudinal  optical  phonons  characterized  by 
their  wave  number  k  and  frequency  to*.  In  terms  of  the  destruction  (creation) 
operators  bk  (bj.)  the  scalar  potential  is  given  by 

<t>  =  (2itfV~'TjMkbk  exp(i k-  r-iukt)  +  h.c.  (1) 

k 

The  coefficients  Mk  arc  at  present  arbitrary.  In  order  to  connect  Mg  to  the  dielec¬ 
tric  properties  of  the  crystal,  we  let  a  charge  travel  through  the  crystal;  i.c.  the 
charge  density  p  is: 

p  =  Qb(r~\)t),  (2) 

and  the  interaction  Hamiltonian  has  the  form: 

//int  -J d/"p0f  CK2rr)3  V~x  Yj  Mkbk  cxp(i&-ut  -  ico*f)  +  h.c.  (3) 

From  time-dependent  perturbation  theory  we  calculate  the  transition  amplitude  to 
the  first  order  in  U„n 

C=(-i//»)  /  dr  <//,„,),  (4) 

—  ©o 

where  <//,,,, >  denotes  the  matrix  element  for  one  phonon  emission.  Standard  cal¬ 
culations  leads  to  the  rate  of  energy  loss  of  the  charged  particle  which  is  given  by: 

-dW/dr  =  E/ico*(2tt)6  V''1{Q2jh')\Mk\'l2nb{o)k  -I t-  u).  (5) 

So  far  our  expression  for  the  energy  loss  is  exact.  Here,  however,  we  intend  to  cal¬ 
culate  A/*,  for  wavenumbers  smaller  than  the  characteristic  lattice  wave  number.  We 
shall  do  this  by  .comparing  a  formula  for  the  energy  loss  obtained  by  using  dielectric 
(or  continuum)  formulation  with  eq.  (5).  We  shall  rewrite  eq.  (5)  considering  only 
long  wavelength  phonons  and  keeping  in  mind  a  cutoff  for  large  wavenumbers.  We’ 
therefore  neglect  the  dispersion  of  the  phonons,  as  usual,  (i.e.  to*  is  replaced  by  coi, 
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the  longitudinal  phonon  frequency)  and  obtain 

-dlV/d/  =  C2it)3  K-*  ftlkQ'WkWlnlVuMun  - k  ■  u),  (6) 

which  is  the  rate  of  energy  loss  of  the  charge  Q  moving  through  the  crystal. 

We  shall  next  take  a  classical  approach  and  calculate  dlt'/df  for  a  charge  Q  trans- 
versing  a  dielectric  medium.  The  charge  density  p  is  given  in  cq.  (2).  The  potential 
induced  by  the  charge  can  be  written  as 

<!>(/•,  r)  =  J  dA  dco  cxp(i/c  •  r  -  itor)  <I>(A,  to),  (7) 

and  similarly  the  charge  density  may  be  written  as: 

p(r.  r)  =  Jtik  dco  exp(i/c  •  r  -  ic or)  p(A,  co).  (8) 

From  the  Poisson  equation  we  have  the  relation  between  potential,  charge  and 
dielectric  function: 

<I>(A,  to)  =  4 rr  p(A,  co)/A2  c(k,  co),  (9) 

where  cqs.  (2)  and  (8)  give: 

p(A,  co)  =  (56(co  -  A  •  o)/8rr3 .  (10) 

Thus,  from  cqs.  (9),  (10)  and  (7),  the  induced  potential  is  now  given  by 

<Kr,  /)  =  J~ dA  du)  cxp(iA  •  r  -  itor)  Q8(co -A- u)/2jr2A2e(A,  to),  ■  (1 1) 

and  the  electric  field  is  given  by  E  =  -V<1>.  Following  Ritchie  [10)  the  work  done 
on  the  charge  is  given  by: 

ir’HZ  «e- « •»■$•/“  (12) 

o 

Here  E0  is  the  induced  field  in  vacuum  which  is  obtained  upon  replacing  c(k,  to)  in 
cq.  (11)  by  unity.  We  have  also  made  use  of  the  reality  condition  e*(A,co)  = 
c(A,  -to). 

We  next  use  the  long-wavelength  limit  of  the  dielectric  function  of  the  crystal 
and  approximate 

c(k,  to)  =»  e(to)  =  €o*  (to2  -  co2)/(co2  -  to2).  (13) 

Merc  to,  is  the  transverse  phonon  frequency,  and  to  has  a  small  positive  imaginary 
constant  to  ensure  that  e(co)  is  the  retarded  dielectric  response.  Using  cqs.  (12)  and 
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(13)  we  obtain: 

dll'  ()2  r  ..  w?  -  w? 

-  —  =  i  / (iA'  j-  ir6(<xj,  -  A:  •  oj. 
d/  rr J  2e„Ai 

Upon  compai  ing  et|s.  (6 )  and  ( 1 4)  we  obtain  [13] 


14 1  (2n)J4nvL  e0)’ 


(14) 


(IS) 


where  e„,  (c0)  is  the  high  (low)  frequency  dielectric  constant  of  the  solid.  This  leads 
to  a  potential 

0(r,  t)  =  K-1,2S(“i)|^-^(i-  -  ~]j  X  l>k  cxp(i k- r  -  iw|/)  +  h.c.,  (16) 

and  to  an  interaction  Hamiltonian  for  the  electron- phonon  system,  given  by 

//,m  =  Jdrpc0=JdreZ)6(r-  rs)<p=  (-i)|^™;,wi (~  ““)] 

Xcxp(iA*  r,  -  iw( /)/;*  +  h.c.,  (17) 

where  rs  represents  the  electronic  coordinates. 

For  the  case  of  semi-infinite  crystal  the  situation  is  somewhat  more  complicated. 
Here  we  lose  our  translational  invariance,  say,  in  the  z  direction.  As  far  as  the  bulk 
properties  are  concerned,  the  semi-infinite  and  the  infinite  crystal  are  identical. 
However  one  should  consider  here  the  reflection  of  the  vibrational  modes  by  the 
sui  face  discontinuity,  which  leads  to  slightly  different  set  of  eigenfunctions. 

We  consider  the  dielectric  approach  for-the  polar  crystal  excitation  spectrum  and 
make  a  Fourier  expansion  of  the  electrostatic  potential 

<l>(r,  t)=J<ik  dw  c.\p(iA  •  r  -  it ot)  <1>(A,  w).  (1 8) 


We  take  the  crystal  for  z  <  0  and  thus  for  z  >  0  the  Laplace  equation  demands  that 

J tlA  dw  cxp(iA  •  r  -  iwf)  k2  <!>(A,  w)  =  0.  (19) 

Hence  AJ<I>(A,  w)  is  an  analytic  function  in  the  upper  Af  plane.  Inserting  this  into 
cq.  (18)  leads  to  the  result 

<l>(r,  /) =  J^dA'i  dw  cxp(iAi-  r  -  iojt)A(klt  to)  c.\p(-AAz),  .  (20) 


forz  >0,  where 

A(klt  w)  =  [(rr/AJ k2  <1>(A,  w)j*I=ikr 
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1  Icrc  /cjl  is  a  two-dimensional  vector  in  the  A--.y  plane.  The  choice  of  the  potential  in 
the  crystal  2  <0  is  somewhat  arbitrary.  However  for  our  purposes  we  choose  the 
potential  to  be  the  sum  of  two  terms  Oiie  of  them  is  a  reflection  of  the  z>0 
potential,  which  makes  <1>  continuous  across  2  =  0.  The  second  term  must  vanish  at 
2=0,  in  order  to  assure  the  continuity  of  <1>  across  the  plane  2=0.  We  thus  write 

00 

<l>(r,  t)=j dAx  du  J  d kz  B(k,  to)  0(— 2)  sin(A.z)  exp(iAx  •  r  -  ico/) 

0 

+  JdAx  dw/lfAx.oj)  cxp(-A.xl2|)cxp(iAx'  r  —  ico/),  (21) 

where  the  integral  on  k.  was  taken  from  0  to  00  to  avoid  double  counting.  Here  cqs. 
(20)  and  (21)  represent  our  general  solution  for  'l»(r,  l),  for  all  space,  in  the  con¬ 
tinuum  approximation.  This  potential  can  be  used  to  describe  electron-phonon 
processes  having  small  momentum  transfer,  i.c.  much  smaller  than  lattice  momen¬ 
tum  for  both  bulk  and  surface  electrons. 

Using  our  choice  of  potential,  one  may  interpret  it  as  having  a  bulk  term  that 
exists  for  2  <  0  and  vanishes  at  the  surface,  and  a  surface  term  which  exponentially 
falls  towards  zero  for  large  values  of  |z|.  Thus,  as  in  the  case  of  metal-insulator— 
scmiconductor  devices  when  electrons  are  trapped  in  the  insulator-semiconductor 
interface,  the  electron- -phonon  interaction  depends  on  the  surface  term  in<l>(r,  t). 
Also  note  that  the  reality  of  <1>  impose  symmetry  relations  on  both  A  and  B.  We 
now  make  contact  with  the  phonon  field  and  express  the  arbitrary  coefficients  A 
and  B  in  terms  of  the  phonon  coordinates.  Here  w*  as  before  is  the  eigenvalue  of 
the  bulk  phonon  having  wavenumber  A .  and  ilkl  is  the  eigenvalue  of  the  localized 
surface  phonon  having  wavenumber  A^.  Our  potential  is  now  given  by: 

<P(r,  t)  =  (2nf  V~x  E  [Mkbk  sin(A.z)  cxp(iAx  •  r  -  i u>kt)  0(-z)  +  h.c.) 

k 

+  (2rr)J/r'  E  lA'x1flx1  cxp(-AxUI)  exp(iAx*  r  -  \9.kyt)  +  h.c.),  (22) 

where  l>k  ( I>1 )  is  the  destruction  (creation)  operator  of  the  bulk  phonon  and  ak]L 
(ttkl)  is  the  destruction  (creation)  operator  of  the  surface  phonon.  The  interaction 
Hamiltonian  is  given  by 

//int  -fi p(a  I)  dr- 

We  arc  however  interested  in  the  surface  modes  and  their  coupling  to  surface  elec¬ 
trons  and  thus  choose  our  p  to  be 

p(r,  t)  =  Q  6(2)  h(rL  -  u x0,  (23) 
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where  rL  and  n,  arc  two-dimensional  vectors  parallel  to  the  surface.  As  before  we 
calculate  the  matrix  clement  (//;,„)  between  no-phonon  and  onc-surface-phonon 
states  and  evaluate  the  energy  loss  rate  to  be: 

-  ~  =QifMi  a,  1 Nh I ' 2 « (ft,  -  kL •  Ujl).  (24) 

In  cq.  (24)  we  have  already  replaced  by  £2(,  the  long  wavelength  surface 
phonon  frequency. 

We  now  calculate,  using  the  classical  approach,  the  energy  loss  of  a  :hargc  mov¬ 
ing  on  the  surface  of  the  crystal.  Here  rather  than  solving  the  potential  in  the  z  <0 
region  for  the  inhomogeneous  problem,  sve  wish  to  consider  an  infinite  crystal. 
1  lowevcr  we  must  impose  an  image  charge  to  account  for  the  dielectric  discontinu¬ 
ity  at  z  =  0.  We  choose  an  induced  surface  charge: 

pf  =  5(z)  o(rlt  t\  (25) 

x  which  together  with  p.  cq.  (23),  and  the  l’oisson  equation  result  in  the  Foutier 
transform  of  the  potential; 


f  Q  ,,  ,  ,  x  o{klt  oij 

Ki’ u) '  F7(TS)ls7,s("  -  ** ' •  ■*> ' *  sri 


where  (j(A‘i.w)  is  the  Fourier  coefficient  of  o(fx,  t).  We  therefore  obtain  for  the 
potential  in  the  2  <  0  region: 

•h(r,  t)  */dA  doj  [^t6(w  -  kj,  ■  vx)  +  —  exp(i*  •  r  -  iut).  (27) 

For  z  >  0  region,  the  potential  is  given  as  before  by: 

•1  >(r,  t)~  f  dAi  dto  /I (Ax,  to)  cxp(-Axz)  c.\p(iAx  *  r  -  icof).  (28) 

We  next  impose  the  usual  boundary  conditions:  (i)  d»  is  continuous  across  z  =  0  and 
(ii)  the  jump  in  D.  is  proportional  to  the  external  surface  charge  density,  and  obtain 
the  following  conditions: 

lAx/e(Ax,  a))j/l(Ax,  to)  =  2 rr  o(Ax,  to)  +  (Q/2a)  5(to  —  Ax  •  i>x), 

C\,klA(kl.  to)  =  -2a  d(AA,  to)  +  (QI2n)  5 (to  -  kL  •  ox).  (29) 

In  cq.  (29),  e(AJt  to)  is  defined  by: 

+  oo 

e(Ax.co)  =  (A3/a)  f  dA.  [A2  c(A,  to)J-’ ,  (30) 


and  cb  is  the  dielectric  constant  of  the  barrier  or  the  oxide.  The  solution  of  cq.  (29) 
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yields 

MkL,  «)  •  J-  »(»-*!•  »J. 

o(*i.  «)  Vl  T  V  --  -v  ■“  8(w  -  kL  •  ut). 
Art2  l  +  e(Ari,  to)  eh 


(31) 


Our  cqs.  (27).  (28)  and  (31)  constitute  the  solution  of  the  electrostatic  field  in  the 
entire  region,  In  order  to  calculate  the  energy-rate  loss  of  the  charge  we  need  the 
field  when  the  crystal  is  replaced  by  its  surrounding  “vacuum”  which  in.  our  case,  is 
the  harrier  having  diclcctiic  constant  ch.  Using  eq.  (30),  c  must  be  replaced  by  cb' 
and  we  obtain  for  2  =  0*  that: 


(E  ■  -  Eq)  •  ux  =  fdAxdto  cxp(iA:1*  r  -  i to/) 

J  \  t  /  Ink  i6b 

Ube(*i,  to)  +  I  J 

Using  the  relation  (e(Arx,  to)]*  =  c(kL,  -  to)  the  energy-loss  rate  is  given  by: 


(32) 


dlf 

dr 


fdk+  f  dco 5  (to  - k±  •  ujj Irnf j). 
J  xkLeb  \c(*lt  to)  +  eb7 


(33) 


We  next  approximate  c(k, to)  by  its  local  (spatially)  approximation  for  longwave- 
length  excitation,  i.e. 


n  ^  _  w5  -  wf 

c(k,  u)se„-; - 

«  -  to, 

Using  eq.  (30)  w:  obtain 

1  toJ  --  to? 
e„  co  -  to? 


and  substituting  it  in  eq.  (33)  results  in  energy-loss  rate: 


(34) 


Upon  comparing  the  quantum  and  classical  expressions  for  the  energy  loss  rates, 
cqs.  (24)  and  (34),  we  conclude  that  the  optic  longitudinal  surface  phonon  fre¬ 
quency  is  given  by 
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ry-l  x  ,-ni'l 
_  Co  +Ct, 

ai~wi  prr,vT  . 

Kw  *  €|»  . 


where  17)  is  always  smaller  than  to}.  We  also  obtain  that: 

ijv  | :  -  _ j  _  _  6i'JL  -  nr\ 

W  (2rr)J  2nAi  +  cft'Xeo'  +£{;')' 

Using  our  e<].  (22)  we  obtain  for  the  potential  on  the  surface  r  0 

a  r  ,-ioV ,  n[  2_«r _  AftjCeJ  -  Co_‘) . T'1 

'  k  U,  C|,(C^'  +  Cf/KCft1  +  c,"‘ )J 

X  exp(~ALU|)  e\p(iAL-  r  -  ii2|/) «fcl  +  h.c.  (37) 

For  electrons  which  arc  bound  to  move  on  the  surface  of  the  crystal  we  have 
P„  =  eT/t>(z)h(r  -  rs). 

X 

where  rs  represent  the  electronic  coordinates.  The  interaction  Hamiltonian  of  this 
two-dimensional  electron  gas  with  surface  phonons  is  given  by 

ft\M  -  f  ll'Pc 
and  reads: 

ii  =  r'1'  ES(-i)f — -  —  tl-l-t€J.rc±l _ 1 1,2 

nm  a  zjIa  feiei)(e-.+£r.Xej1+e-.)J 

X  espiiAi*  rs  -  \Sl{t)akl  +  h.c.  (38) 

The  interaction  Hamiltonian  as  given  by  etj.  (38)  repicsents  the  inteiaction  of  a 
two-dimensional  electron  gas  with  optic  longitudinal  vibration  on  the  polar  crystal 
interlace. 


3.  Conclusion  and  discussion 

Using  the  dielectric  formulation  we  were  able  to  detenninc  the  interaction 
between  electrons  bound  to  move  on  a  surface  of  a  polar  semiconductor,  and  lattice 
vibrations.  This  interaction  is  the  result  of  electrons  absorbing  or  emitting  surface 
phonons  (optic-longitudinal),  '(he  interaction  potential  in  an  electron- phonon  col¬ 
lision  is  given  by: 


v  =  ..j  .  /,sMc~  "  c«')  __T2 

K>  U’l  C|,(C'  <  e,~' )  (co*  + 
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The  surface  phonon  frequency  was  also  obtained  and  is  defined  in  cq.  (35)  in  terms 
of  tlte  bulk  longitudi  ir.i  phonon  frequency.  The  electron-lattice  time-independent 
interaction  Hamiltonian  can  be  written  as 

Hint  =  '1  ‘ 1,2  S  vkl  pklakl  +  h.c.,  (40) 

where 

Pfcj.  =  Z)cxp(iAcj.-  rs) 
s 

is  the  l-ouricr  transform  of  the  electronic  density  operator  which  can  be  repre¬ 
sented  by 

S  Cp+kicp , 
p 

where  cp  (cj,)  arc  the  usual  destruction  (creation)  operators  for  the  electron  gas. 
The  interaction  Hamiltonian,  cq.  (40),  is  the  appropriate  one  to  use  when  electron- 
phonon  processes  are  considered,  for  electrons  that  arc  bound  to  move  on  an  insu¬ 
lator-semiconductor  (polar)  interface. 
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We  calculate  the  high-frequency  conductivity  lor  a  two-dimensional  electron  gas  interacting 
with  optical  lattice  vibrations.  The  collision  trcqucncy  and  the  optical  mass  normalization  are 
obtained.  The  frequency  dependence  of  the  (inverse)  collision  time  is  presented. 


1.  INTRODUCTION 

When  a  device  of  metal-insulator-semiconductor 
structure  is  biased,  a  large  accumulation  of  electrons 
or  inversion  layer  is  produced  on  the  surface  of  the 
semiconductor.  When  the  energy  bands  are  bent,  as 
in  the  case  of  strong  inversion,  a  potential  welt  ap¬ 
pears  on  the  insulator-semiconductor  interlace  which 
localizes  the  motion  of  the  electrons  normally  lo  the 
surface.  On  the  other  hand,  the  electrons  are  free  to 
move  in  the  plane  of  the  surface.  Under  strong 
enough  bias,  the  localization  of  the  electrons  in  the 
direction  normal  to  the  surface,  which  is  character¬ 
ized  by  the  spread  of  their  wave  function,  is  much 
smaller  than  their  average  distance  within  the  plane. 
Wc  are  thus  permitted  to  use  a  two-dimensional 
model  for  the  electron  gas  in  which  the  electrons  arc 
free  to  move  only  in  the  plane  of  the  semiconductor 
surface.1 

In  treating  the  response  of  the  electron  gas  in  our 
system  lo  a  dc  electric  field  (transport  theory)  or  an 
ac  electric  field  (optical  properties),  collision 
processes  of  electrons  with  tons,  impurities,  or  acous¬ 
tic  phonons  have  been  mainly  considered.2  These 
can  explain  results  obtained  in  semiconductors  such 
as  Ge  and  Si.  However  when  compound  semicon¬ 
ductors  are  used,  such  as  InSb,  GaAr,  etc.,  which  are 
partially  ionic,  the  interaction  of  the  electrons  with 
longitudinal-optic  phonons  becomes  important  and 
cannot  be  neglected.-'  A  similar  situation  exists  for 
bulk  electrons  in  compound  semiconductors.  The  ef¬ 
fect  is  particularly  notable  at  low  temperature  for  opt¬ 
ical  absorption  when  the  photon  frequency  w  exceeds 
the  longitudmal-optic-phonon  frequency  U(.  Here  a 
new  channel  of  absorption,  i.e.,  a  final  state  with  an 
excited  phonon,  is  operative  and  influences  strongly 
the  optical-absorption  coefficient.  In  Sec.  II  we 
present  the  theory  for  the  high-frequency  conductivi¬ 
ty  of  a  two-dimensional  electron  gas  with  polar  optical 
phonons. 

In  treating  the  electron-phonon  interaction  one 
should  note  that  the  electrons  moving  on  the  surface 


of  the  semiconductor  will  interact  most  strongly  with 
surface  rather  than  bulk  phonons.4  Surface  phonons 
produce  a  large  electrostatic  potential  in  the 
insulator-semiconductor  interface.  The  potential  de¬ 
creases  exponentially  to  zero  away  from  the  interface. 
The  etgenfrcqucncies  of  the  surface  phonons  and 
their  interaction  with  the  electrons  depend  on  bulk 
properties  as  well  as  on  the  dielectric  properties  of 
the  insulating  barrier. 


11.  CALCULATIONS  OF  THE  CONDUCTIVITY 

Our  electron-phonon  system  is  described  by  the 
Hamiltonian 

//-//«+//,  .  .  U) 

where 

ll0  **  £  b,  (2) 

p  i 

and 

//l“-~  viap*qap'-tap'ap 
Ln  p.p'.i 

+  ~7uT  +ll,c.)  .  0) 

"  pi 

Here  tp  is  the  kinetic  energy  of  an  electron  having 
momentum  p.  to,  is  the  wavc-number-depcndcnt 
longitudinal-optic  frequency  and  a }  (/>,,  />„’)  arc 

the  destruction  and  creation  operators,  respectively, 
of  the  electrons  (phonons).  The  coupling  term 
v,  ~2ne!lti  is  the  Fourier  transform  of  the  Coulomb 
interaction  for  planar  electrons  and 

Cf— fa.,  [«»(*;'  +€»'),fl(«0 '+«*', '"I-'/1 

</ 

represents  the  electron-longiludinal-optic-phonon  in- 
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teraction. 

To  evaluate  the  conductivity  we  start  from  the' 
Kubo  formula  which  reads 

("(w)  -  J*o  </X  (j  (r  -  (X) -j  (0))  ,  (4) 


where  to  is  the  frequency  of  the  applied  field,  J(0)  is 
the  Fourier  transform  of  the  current  operator  for 
wave  number  equal  to  zero,  and 

J(T)  1(0)  <*-•"'  .  (S 

The  statistical  average  of  an  operator  0  is  given  by 

<0>-Tr(eflUU"w-w,0)  .  (6 

where  //  is  the  total  Hamiltonian  of  the  system  and 
(1  is  defined  by 

e-fln-Tr(c^'-">)  .  .(7 

Here  /a  and  N  are  the  chemical  potential  and  the 
number  operator,  respectively,  and  fi  is  the  inverse 
of  the  temperature  in  energy  units. 

In  order  to  render  F.q.  (4)  in  a  more  convenient 
form,  we  integrate  by  parts  and  obtain 

IT  -  <T0  +  IT,  .  (8 


In  Eq.  (10)  the  square  brackets  denote  the  commuta¬ 
tor. 

Calculations  of  the  current-current  correlations  and 
the  conductivity  in  the  three-dimensional  case  have 
been  worked  out  in  much  detail  and  arc  well  docu¬ 
mented.5  The  calculations  for  the  two-dimensional  sit¬ 
uation  is  remarkably  similar  and  we  shall  present 
only  the  final  results.  We  evaluate  the  conductivity 
treating  electron-phonon  collision  within  the  Dorn  ap¬ 
proximation  (high-frequency  conductivity),  however 
treating  the  self-consistent  field  of  the  fluctuating 
electron  gas  exactly  in  the  random-phase  approxima¬ 
tion  (RPA).  Our  expression  includes  the  full  dynam¬ 
ic  screening  of  the  electrons.  We  however  will  con¬ 
sider  here  a  weak  electron-phonon  interaction  and  ig¬ 
nore  completely  the  renormalization  of  the  phonon 
spectrum  and  line  broadening  by  the  electron  density 
fluctuations.  Our  result  for  the  conductivity  reads 


+  -  3  J* J  dx coth  — r-^— 

2wV(27t)1  j  '  4w  2  «,(*+« 


x  — —  ((),(*  +«j)  -(?,(x)IUVx  +<a)  -£>,(*)) 
^  ) 


- 7p-[(?,U+<a)-(?/U)l|D,(x+cu)-D/U))  .  (11) 


Here  the  dielectric  function  is  given  by 
«,(*)  - 1  -  v„Q,(x)  , 

where  Q,  the  density  fluctuation  is  defined  by 

<?,<*>  -TTw/^P - — -  • 

(27f)'  J  6,+,/J -ip-'H- X  -IT) 


Dp(x)  represents  the  phonon  propagator  and  is  given 

by 

DAx)~-p*-rr  .  (14) 

•  x‘  -  <D, 

To  render  our  result  for  the  conductivity  in  a  more 
transparent  form  we  shall  identify  the  effect  of  the 
electron-phonon  collision  with  the  Drude  form  of  the 


conductivity.  We  write  for  <r, 


(/»  +8/h  )(w  +  /v) 


where  hm,  the  mass  renormalization,  and  v,  the  colli¬ 
sion  frequency,  are  small  quantities.  We  therefore 
identify  them  as 

M-(8jrJw/H/i)-1  Re  J*  r/qi?J|Q|J-^— 

x  P  f  **  dx  coth  F  .  (16) 

J  2 

Sin-- Im  f  r/qt/2|C,P~ - 
J  4ff 

xP  f+"  t/xcoth  F  .  (17) 
»  2 
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«,(x+oj)  t,(x) 


[@,(x  +oj)  -  @?Or))[D,(x  +«)  -Z),(x)) 


-—^tC),(x+a.)-<?,*0.)l[D,(x+a,)-Df*U)l|  .  (18) 

In  order  to  calculate  the  collision  frequency  i>  we  shall  use  two  simplified  assumptions  which  arc  well  justified. 
First  we  omit  all  phonon  lifetime  effects,  thus  ignoring  the  narrow  spread  of  the  phonon  spectrum  for  any  save 
number  q.  Second,  since  the  momentum  transfer  q  is  of  the  order  of  the  Fermi  momentum  which  in  turn  is 
much  smaller  than  lattice  momentum,  we  ignore  entirely  the  dispersion  of  the  phonons  and  replace  w,  by  If/,  the 
longitudinal  vibration  frequency.  After  some  algebra  one  obtains 

v  -  (8rrumm)-1  f  dq  q2 

x  colh  +coih^”(ta  —  II/)  lm — - — —  +  coth  +coth^(w  +  II/)  Ini — ~~r  .  (19) 

2  2  f,((u-Il /)  2  2  «,((u+li /) 

For  most  practical  situations,  the  temperature  is  smaller  than  the  Fermi  energy.  We  therefore  are  justified  in  tak¬ 
ing  the  dielectric  function  at  zero  temperature.  The  statistical  factor  can  also  be  taken  at  the  zero-temperature 
limit  for  kaT  <  All/.  We  thus  obtain  for  the  collision  frequency  at  T  =  0, 


i>  -  (Atu>iiiw)~'  f  dq  q3  ■  ■-  lm  — ; — ■  ■■■  7  f)(w  — 11,)  . 

J  V,  «,(tO-ll/) 


Here,  as  expected,  a  threshold  for  absorption  occurs  at  i«*=  Il(.  The  collision  frequency  rises  for  w  >  11/  anil,  as 
we  shall  see  later,  behaves  as  <o'il2  at  high  frequencies.  1  he  effect  of  finite  temperature  when  A „  T  <  AH/  would 
be  mainly  to  round  off  the  sharp  threshold  at  o»«*  II/.  Numerical  integration  of  Eq.  (20)  will  be  presenteu  later. 
The  mass  renormalization  8m  can  bo  cast  into  a  simpler  form  after  some  algebra  and  reads 


8m  -  (8  TTJa>Jrt)~*  f  tlq  q2  f  dx  coth 

J  [  u,  l  2 


x  lm — 77— r-  Rc(Z2,(x  +o>)  +  D,[x  -w)  -2D, (*)) 

^  \X  / 


+  im D*(x)  Re 


(,{x+Ol)  f,(x  —  o>)  t,(x) 


Numerical  calculations  of  8m  as  a  function  of  fre¬ 
quency  arc  more  difficult  than  for  v.  Here  the  vari¬ 
ous  terms  in  the  integrand  have  alternating  signs, 
thus  large  cancellations  do  occur.  Let  us  determine 
the  qualitative  behavior  of  8m  as  a  function  o>.  Merc 
for  u— ‘oo,  t,(x  ±01)  may  be  replaced  by  unity  and 
D,(x  ±  01)  asymptotically  approaches  zero  as  o'2.  We 
thus  conclude  that  for  w  — »  the  integral  in  F.q.  (21) 
approaches  a  constant  and  hence  8m  —  o>'2.  On  the 
other  hand  for  ai— 0  (below  the  longitudinal  phonon 
frequency)  the  integrand  in  Hq.  (21)  behaves  as  o>2, 
and  thus  6m  becomes  a  constant  when  o>  approaches 
zero.  A  detailed  behavior  of  8m  as  a  function  of  o> 
will  not  be  presented  here. 

We  go  back  now  to  our  Eq.  (20)  and  evaluate  the 
collision  frequency  i>.  Here  we  consider  the  situation 


where  the  Fermi  energy  is  larger  than  kaT and  there¬ 
fore  arc  justified  in  using  for  c,(<»-  11/)  the  zero- 
temperature  dielectric  function.  Our  result  reads 

—  -  {(Ac)-jr  f"  dz  z~  <)(  11  -  ill)  , 

o>i  1  11  tU  11  - 11,) 


where  01,  is  the  longitudinal-optic-phonon  frequency, 
11,  is  the  surface  phonon  frequency, 

A  *»  [lc2in/ K2(ttii),,2\  is  the  plasma  parameter  of  the 
two-dimensional  electron  gas  and 
{ -1  _  f  -1 

is  the  dielectric  form  factor  for  the  interaction. of  the 
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electrons  with  surface  phonons.  Here  0  -a»/4«/-  and 
z  -q/2kr  arc,  respectively,  the  frequency  and  wave 
numbers  normalized  with  respect  to  the  Fermi  energy 
and  momentum.  Also  11/  is  the  surface  phonon  fre¬ 
quency  similarly  normalized.  The  integrand  in  Eq.  (22) 
is  given  by 


F(z.x)-D .  1  -  z— - 
z 


0.  if|z+-|ssi 

V 

1.  if  |;  ±  -j|  <  1 


(23) 


and 

t(z,x)  -  1  +  —  G (z, x)  +  /-q~F(z,x) 
z  2  z 

where 


G(z.x) 


0.  if|z±£|*SI  , 

■  :r  l,  +  £|  - 


(24) 


F.quation  (22)  can  now  be  integrated  using  InSb 
parameters:  «j,  =  24.4  mcV,  «0°  17.9,  15.7,  and 

as  an  example  for  the  oxide  layer,  tb  -  10  is  taken  as 
a  reasonable  value.  The  Fermi  energy  in  our  calcula¬ 
tion  is  taken  to  be  100  mcV.  In  this  case  we  obtain 
for  the  surface  phonon 

,.-i  '« 

Sh-OJ,  - b-r  “0.95o»|  . 

+  <» 


which  is  almost  equal  to  <u,.  Our  result  for  v  as  a 
function  of  the  normalized  frequency  is  presented  in 
Fig.  I.  The  collision  frequency  v  is  zero  for  frequen- 


FIG.  I.  Normalized  collision  frequency  P-2i7 xCw,  as  a 
function  of  a  -  (w  -  0,)/4w/-  for  InSb  at  zero  temperature, 
for  X  -0.4  and  X  -  I. 


cics  below  the  surface  phonon  frequency.  When  the 
light  frequency  exceeds  the  surface  phonon  frequen¬ 
cy,  v  increases  fast  up  to  a  constant  value  and 
remains  so  until  w>»4a>,'.  Thereafter  v  approaches 
asymptotically  to  a  constant  times  o>~sl\  for  large  w. 
This  asymptotic  form  of  v  is  obtained  when  screening 
effects  are  omitted,  i.e.,  when  in  Eq.  (22)  we  replace 
<(2,  fl  -  12/)  by  unity. 

In  conclusion,  we  have  calculated  the  collision  fre¬ 
quency  for  a  two-dimensional  electron  gas  interacting 
with  optical  lattice  vibrations.  Our  result  for  the  col¬ 
lision  frequency  as  a  function  of  the  photon  energy  is 
plotted  in  Fig.  I.  We  find  that  metal-oxide- 
semiconductor  devices,  made  up  of  polar  semicon¬ 
ductors,  will  exhibit  a  jump  in  the  collision  frequency 
when  the  photon  frequency  exceeds  the  optical-pho¬ 
non  frequency. 
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High-frequency  conductivity  of  a  two-dimensional,  two-component  electron  gas 
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A  first-principles  calculation  of  the  high-frequency  conductivity  for  a  two-dimensional,  two-component 
electron  gas  is  presented.  Both  frequency  and  temperature  dependence  arc  considered.  The  influence  of 
electron-electron  collisions  on  the  two-dimensional  plasmon  is  discussed  in  relation  to  possible  experiments 
in  MOS  structures. 
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I.  INTRODUCTION 

Recently,  there  has  been  considerable  experi¬ 
mental  and  theoretical  interest  in  the  physics  of 
electrons  in  so-called  metal-oxidc-semiconductor 
(MOS)  structures.1  In  these  devices,  electrons 
are  bound  by  the  application  of  an  applied  elec-  ' 
trie  field  to  the  oxide-semiconductor  interface. 
The  potential  they  experience  depends  in  detail 
on  Ihe  field,  the  carrier  concentration,  and  on 
the  dielectric  properties  of  the  sandwich,  it 
suffices  to  say  that  the  effect  of  this  potential  is 
to  quantize  motion  perpendicular  to  the  struc¬ 
ture  into  a  set  of  discrete  levels  or  subbanos  with 
spacing  of  the  order  of  100  "K.  The  motion  in  the 
plane  of  the  interface  is  nearly  free-electron-like 
with  an  effective  mass  characteristic  of  the  three- 
dimensional  (3-D)  band  structure.  At  low  tem¬ 
peratures,  this  sytem  behaves  in  many  ways  like 
an  almost  ideal  2-D  electron  gas.  The  fact  that 
its  density  may  be  varied  over  many  orders  of 
magnitude  by  simply  changing  the  applied  elec¬ 
tronic  field  has  made  it  an  interesting  system  for 
the  study  of  Coulomb-interaction  effects  among 
electrons. 

The  high-frequency  conductivity  of  such  a  sys¬ 
tem  reflects  many  of  its  important  physical  prop¬ 
erties.  In  fact,  it  lias  been  known  for  some  time 
that  the  cyclotron  resonance  line,  or  lires,  in 
such  systems  depends  in  a  rather  complicated  way 
on  the  carrier  density,  the  frequency  of  the  ex¬ 
periment,  the  temperature,  and  the  orientation  of 
the  interface  relative  to  the  crystal  axis.2  Several 
attempts  to  understand  the  behavior  of  such  cyclo¬ 
tron-resonance  experiments  within  the  framework 
of  an  interacting  2-D  electron-gas  picture  have 
been  made.5  All  such  discussions  except  one'1  are 
based  on  a  picture  wherein  a  single  group  of  car¬ 
riers  interacts  with  one  another  and  with  impur¬ 
ities.  These  calculations,  roughly  speaking,  give 
an  adequate  description  of  the  cyclotron  resor  nee 


linewidth,  the  presence  of  subharmonics,  and  in 
some  cases,  the  position  of  the  line  as  a  function 
of  carrier  concentration. 

One  set  of  experiments  on  Si(100)  MOS  struc¬ 
tures  have  shown  that  the  temperature  dependence 
for  5  T<  <  T<  80  K  is  anomalous.  If  we  take  this 
one  experiment  at  faee  value,  we  are  forced  to 
conclude  that  the  single  resonance  line  moves 
from  a  position  determined  by  an  effective  mass 
of  0.2  mc  to  an  effective  mass  near  0.35m,  as  the 
temperature  is  raised.  Attempts  to  understand 
this  effect  within  the  conventional  many -body  pic¬ 
ture  discussed  in  the  preceding  paragraph  have 
failed  miserably.  In  fact,  estimates  of  the  tem¬ 
perature  effect  which  rely  on  a  single-species 
interacting  Coulomb  liquid  suggest  that  the  mass 
should  become  lighter  as  the  temperature  is  . 
raised.5 

Fortunately,  for  Si(100),  the  situation  is  some¬ 
what  more  complicated  than  our  simple  discussion 
would  indicate.  In  this  case,  the  multiellipsoidal 
nature  of  the  3-D  band  structure  gives  rise  to  two 
sets  of  subbands  in  the  2-D  case.  One  of  these, 
the  lowest,  had  a  2-D  mass  of  roughly  0.2w,  while 
the  second  set  has  a  cyclotron  mass  of  roughly 
0.4m,.  Appel  and  Overhauser*  have  recently  used 
this  fact  to  discuss  the  temperature  dependence 
of  the  cyclotron  resonance  line.  Using  a  set  of 
mncroscopic  transport  equations  for  the  current 
along  with  an  effective  electron-electron  relaxa¬ 
tion  time,  they  were  able  to  show  that  one  could 
qualitatively  understand  the  CR  experiments  as 
a  collisionally  averaged  line  for  the  two  groups 
of  carriers.  Of  course,  there  is  still  some  ques¬ 
tion  of  whether  this  is  the  correct  interpretation 
of  the  data.  The  idea,  however,  is  an  interesting 
one  and  seems  to  deserve  further  microscopic  an¬ 
alysis. 

In  this  paper,  we  present  a  first -principles 
calculation  of  the  high-frequency  conductivity 
a(w)  of  a  two-component  degenerate  system.0  The 
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calculation  is  strictly  valid  to  lowest  order  in 
\/2  =  r,*chn/ti2kF,  where  »i  and  kF  are  the  mass 
and  Fermi  wave  vector  of  the  light  carrier.  It 
properly  includes  all  the  dynamic  screening  ef¬ 
fects,  The  results  will  be  useful  for  interpreting 
the  far-infrared  Drude-typo  conductivity  of  the 
inversion  layers.  As  we  shall  see,  correlation 
effects  which  appear  in  the  conductivity  will  in¬ 
fluence  the  position  and  linewidth  of  the  2-D  plas- 
mon  resonances  which  have  been  observed.1 


II.  PORMULATION 

The  Hamiltonian  of  our  system  is  given  by 

//=  23  o~  l’},a,%a  +  %  23,  vlaUZ,*aVZ,a,nV,«‘ni,°!- 

p  a  J  J  o 

(D 

In  Eq.  (1),  a  (n|  a)  is  the  destruction  (creation) 
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operator  of  electrons  of  type  "a"  having  a  momen¬ 
tum  p  and  vt=2nc2/<i  is  the  Fourier  transform  of 
the  Coulomb  interaction  between  electrons  con¬ 
fined  to  move  in  a  plane. 

Rigorous  solution  for  the  conductivity  for  the 
full  range  of  frequencies  is  difficult  to  obtain. 
However,  for  the  high-frequency  conductivity, 
where  wx>  1  (t  some  effective  collision  time),  a 
solution  for  o(w)  was  obtained  in  the  3-D  case.8 
Here,  collisions  between  particles  having  large 
momentum  transfer  are  treated  within  the  Born 
approximation.  However,  interaction  effects 
within  the  Hartree  approximation  are  considered 
to  all  orders.  We  thus  treat  the  local  field  effect 
exactly,  i.e.,  in  the  full  random-phase  approxima¬ 
tion  (RPA). 

The  calculation  of  the  conductivity  in  the  2-D 
case  can  be  carried  out  in  a  similar  fashion  to  the 
3-D  case,  which  is  well  documented.  The  final 
result  of  such  a  calculation  is 


(2) 


The  correction  to  a  [second  term  in  Eq.  (2)]  vanishes  when  ni=M.  This  is  a  consequence  of  the  general 
theorem  that  the  long-wavelength  conductivity  is  the  noninteracting  conductivity  when  the  particles  all 
have  the  same  mass,  i.e.,  conservation  of  momentum  implies  that  the  total  current  is  conserved.8 

In  Eq.  (2),  »,»»»,  (N,M)  are  respectively  the  "light"  electron  ("heavy  electron)  density  and  mass,  and. 

11  —  u)/AEF  where  EF  is  the  light-electron  Fermi  energy.  Also,  z  =  q/2kF  represents  the  momentum  trans¬ 
fer  in  the  collision  process.  The  retardation  effects  are  represented  by  the  integration  over  x  (P  is  the 
principal  value),  which  is  the  Fourier  transform  over  the  time  in  our  normalized  energy  units. 

The  function  /•’  is  quite  complicated  and  has  a  statistical  factor  depending  on  ,1  =  4 EF/kBT  where  kB  is  the 
Boltzmann  constant  and  T  the  temperature.  The  rest  of  the  terms  in  our  expression  for  F  represent  the 
long-wavelength  current  fluctuations  in  terms  of  the  convolution  of  the  density  fluctuations  of  the  light  and 
heavy  electrons.  It  reads 


F-i coth  €;1(x  +  u>){«;,(.v)[G,(x+ w)-G,(x)l[/t(v  + w)-Jt(x)] 


-t;-'(.v)[Gf(x+w)  -G/(x )][</,(* +  w)  -J»1} . 


(3) 


Here  G,{x)  [</,(*)]  is  proportional  to  the  density 
fluctuations  of  the  light  (heavy)  electrons  and  is 
given  by 


Gt(x)=G*(x)+iG'(x) , 
where 

c;<*)=i-£[(*-*/z),-i],/i 

-^«z  +  x/z)*- If'*, 
0,  if  | z  ±x/z  |  «  1 


(4) 


(5) 


-63- 


and 


<?:(*>=£  K(i -<*-*/«  )T2 

-D,[l-(z+x/z)*]*/*|  t 
Dt(0,  if  |z±x/z  |*  1 
l  1,  if  jz±x/z  |<1.. 


Similarly, 

</,(*)  =  W+ »'«/'(*) 
and 


(6) 

(7) 

(8) 
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The  dielectric  function  e,(v)  is  defined  for  Our 
system  as 

+  .  (9) 

Our  result  for  the  high-frequency  conductivity, 
Eq.  (2),  to  lowest  order  in  the  plasma  parameter 
is  given  in  terms  of  F,  Eq.  (3),  and  represents 
the  interaction  between  the  density  fluctuations  of 
the  light  and  heavy  electrons  including  their 
dynamical  screening.  Although,  at  first  sight,  our 
result  looks  complicated,  it  represents  a  simple 
pysical  process  for  absorption.  Since  G  and  J  are 
proportional  lo  (he  iight  and  heavy  electron  density 
fluctuation,  respectively  fsee  Eq.  (9)1,  the  final 
state  after  absorption  includes  a  light  electron- 
hole  pair  as  well  as  a  heavy  electron-hole  pair 
The  absorption  process  is  as  follows:  The  incom¬ 
ing  photon  excites,  for  example,  the  light  electron- 
hole  pair  which,  in  turn,  via  the  Coulomb  interac¬ 
tion,  excites  a  heavy  electron-hole  pair.  Similar¬ 
ly,  the  photon  may  excite  a  heavy  electron-hole 
pair  which  is  followed  by  the  excitation  of  a  light 
electron-hole  pair  via  the  Coulomb  interaction. 

Note  that  there  is  no  final  state  with  either  two 
light  electron-hole  pairs  or  two  heavy  electron- 
hole  pairs.  These  final  states  do  not  contribute  to 
the  absorptive  part  of  the  long-wavelength  conduc¬ 
tivity.  Also  note  that  in  the  limit  of  »i  =.U  (where 
no  distinction  between  light  and  heavy  electron  oc¬ 
curs)  there  is  no  contribution  to  the  absorption  due 
to  correlation  effects.  Thus  our  first-principles 
theory  accounts  correctly  for  the  appropriate  final 
state  and  also  takes  into  account  the  complicated 
dynamical  screening  represented  here  by  the  two 
factors  of  <•',  which  makes  ICq,  (3)  look  complicat¬ 
ed. 

The  exact  evaluation  of  the  double  integral  in 
Eq.  (2)  is  quite  complicated.  We  may  use  a  sim¬ 
plified  form  for  F  which  proved  to  be  numerically 
acceptable  in  the  three-dimensional  case.  We 
found  that  when  the  frequency  w  was  smaller  than 
w^,  some  mean  plasma  frequency  dynamical 
screening  effects  arc  not  very  important.  In  this 
case,  it  was  possible  to  replace  e4(.v)  by  its  static 
value.  In  this  limit,  Eq.  (3)  becomes 

F = {K('v + «•>)  “  <;«('  )1  K(-v  +  w) -./,<v)l 

-  (G,(*  +  u>)  -G4(.v)][J,(x  +  w)  -J;(x)]} . 

(10) 

There  is  still  an  interesting  physical  question 
present  in  our  choice  of  «4(0).  Formally  [see  Eqs. 
(5)  and  (8)j, 

€4(0)=1  +  j(1  +  .W/»i).  (11) 


We  have  static  screening  from  both  carriers  and 
that  screening  is  independent  of  the  number  of 
those  carriers.  This  is  true  because  the  static 
screening  length  in  a  2-D  degenerate  system  is 
independent  of  the  number  of  carriers.  We  will  be 
interested  in  the  case  where  there  is  a  variable 
number  of  heavy  carriers.  In  particular,  the  ap¬ 
propriate  response  time  (plasma  frequency)  for 
the  small  number  of  heavy  carriers  will  be  small 
compared  to  the  applied  frequency.  In  this  case, 
it  is  a  much  better  approximation  to  neglect  the 
screening  from  the  heavier  carriers  altogether, 
i.o.,  we  pick  e4(0)2  1  f  (X/a)Gj(0). 

After  some  algebra  and  use  of  the  symmetry 
properties  of  GR>r  and  JR'‘  and  the  fact  that  we 
integrate  over.v,  we  find  that  the  reactive  or  real 
part  of  F  can  be  written  as 

F*  =  coth(y){G[(.v)[jJ(x  +  Sl)  +  J,*(x  - 1 <)  -  2J»1 

+  (G~J)}  (12) 

(in  the  second  term,  G  and  J  are  interchanged). 
The  imaginary  or  absorptive  term  can  be  cast  into 
the  form 

Fr~  |coth!j-coth*-(.v  -iI)]g'(SI  -x)J‘t(x).(  13) 

Using  Eqs.  (12)  and  (13),  the  conductivity  Eq. 

(2)  can  be  written  as 


where  the  complex  frequency-dependent  collision 
frequency 


We  have  evaluated  Eq.  (15)  numerically  for  a 
range  of  frequency  and  temperatures  for  vi/M 
=  1/2,  n/N-  4,  and  A=l.  We  find  that  the  real 
part  of  the  collision  frequency  (F{)  can  be  approxi¬ 
mated  by 

t»sr„(ll)  =  2.4^«'+40^yj.  (16) 

The  functional  form  of  the  results  is  to  be  expect¬ 
ed.  The  ST*  and  T1  dependence  is  simply  related 
to  the  smallness  of  the, phase  space  available  to 
pairs  of  colliding  particles  having  an  energy 
transfer  11.  Our  result  for  the  collision  frequency 
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was  obtained  using  microscopic  calculations  for 
the  response  of  the  plasma  to  the  radiation ‘field. 

It  incorporates,  in  principles,  the  dynamical  cor¬ 
relation  effects  of  the  self-consistent  field  gener¬ 
ated  by  the  charged  particles.  However,  for  small 
frequencies  and  temperatures  when  only  low-fre¬ 
quency  density  fluctuations  are  important,  our 
numerical  results  are  essentially  similar  to  those 
presented  by  Appel  and  Overhauser.  The  numer¬ 
ical  coefficients,  of  course,  depend  on  all  the 
details,  such  as  the  number  of  particles,  mass 
ratios,- etc.  It  seems  inappropriate  to  study  the 
dependence  on  these  parameters  without  some 
experimental  data  with  which  to  compare  them. 

The  imaginary  or  reactive  part  of  the  collision 
frequency  is  quite  different.  It  is  not  very  temper¬ 
ature  dependent.  It  has  a  magnitude  at  small  Cl 
which  is  comparable  to  uR  at  11=0.01.  Thus,  for 
any  physically  interesting  situation,  it  will  be 
quantitatively  unimportant.  We  will  neglect  it. 
Thus,  our  first -principles  calculation  of  the  con¬ 
ductivity  shows  us  that  for  frequencies  high  com¬ 
pared  to  the  collision  frequency  itself,  but  low 
compared  to  the  Fermi  energy,  one  must  use  a 
frequency -dependent  real  v  [see  Eq.  (14)].  At 
lower  frequencies  (low  compared  to  v  itself),  a 
macroscopic  semiphenomenological  treatment  of 
a  (see  Appendix)  similar  to  that  used  in  Kef.  1 
leads  to  a  simple  form  for  a,  for  all  frequencies, 
i.e., 

A  -  -\f]  W />»)(!  -iii/M?  v  \ 

a  1  mw  V+Af  till  + 1,  w  N\  / n  M\iCl-v  r 

\  V+AW\N+m)  / 

(17) 


For  low  enough  frequencies,  the  conductivity  de¬ 
velops  a  reactive  part  as  a  result  of  correlations, 
This  reactive  part  will  lead  to  a  shift  of  the  ob¬ 
served  2-D  plasma  frequency.  For  values  of  the 
parameters  used  in  the  numerical  evaluation  and 
assuming  Hu)<kBT<EF, 


~  •  I,  x  1  f  22.5fa  T/ErY  ]) 
r~  ,0#(1  27  U— 22.5(*8r/EF# 


(18) 


There  will  be  about  a  4%  shift  in  the  plasmon  fre¬ 
quency  for  u<{22.5)(kBT)'/EF anda  Lorcntzian- 
type  contribution  to  the  linewidth  as  a  function  of 
T  or  w  of  the  same  order  of  magnitude.  Doth  such 
effects  should  be  observable,  . 
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APPENDIX:  CALCULATIONS  OF  CONDUCTIVITY  FOR 
CONSTANT  COLLISION  TIME 


We  consider  here  an  interacting  system  of  light 
and  heavy  electrons.  As  is  well  known,  a  long 
wavelength  radiation  field  cannot  be  absorbed  due 
to  the  collision  between  electrons  having  the  3ame 
mass.  Thus,  no  relaxation  process  involving 
either  of  the  species  is  possible.  The  collision 
between  the  light  and  heavy. electrons  is  effective 
as  long  as  their  relative  momentum  is  not  zero. 
For  equal  numbers  of  electrons,  the  relative  mo¬ 
mentum  is  given  by  (A//>  -  mP)/(m  +A/)  where  low¬ 
er  case  and  capital  letters  refer  to  the  different 
densities  n  and  ,V  of  the  two  types  of  electrons. 
The  kinetic  equation  reads 


)  itNwAf  /  p  P  \  1 
ut '  mu + NAf  \  mu  NM/r  ’ 


(Al) 


The  second  term  on  the  right-hand  side  of  Eq. 
(Al)  represents  the  relaxation  of  relative  momen¬ 
tum  for  a  constant  electron  collision  time.  Sim¬ 
ilarly,  we  obtain 


P  = 


iiNiiiM  /  P  p  \1 
uni+N M  \  NM  mu  )t  ' 


(A2) 


The  solution  for  the  velocities  v-p/m  and  V-P/M 
for  a  field  of  the  form  Ee,ut  involves  some  al¬ 
gebra  but  is  not  difficult  to  obtain.  Given  the  ve¬ 
locities,  we  can  evaluate  the  current  density 
j  -  c{v+  V)  and  use  a=;/£  to  obtain  our  final  re¬ 
sult.  The  conductivity  (l/r  =  v)  is 


x 


(A/ /»<)(!  -iii/M)1  v  \ 
/,  in  N\  / n  AAfn-pf 

\Urtn)(rm)  ) 


(A3) 


In  other  words,  Eq.  (17)  appropriately  extends 
our  high-frequency  result  to  the  collision-dom¬ 
inated  regime,  i.e.,  it  is  consistent  with  kinetic 
theory  when  a  constant  collision  time  is  consider¬ 
ed. 
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